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ABSTRACT 


The  work  presented  in  this  report  falls  into  two  classes, 
namely  on  the  reflectivity  of  solar  infrared  from  high  altitude  clouds, 
and  on  the  effect  of  the  energy  reflected  from  these  clouds  on  the 
statistical  performance  criteria  of  a  satellite-borne  missile  detection 
system.  The  scattering  of  infrared  from  clouds  of  spherical  water  drop¬ 
lets  has  been  treated  by  Mie  theory.  For  wavelengths  in  the  range 
2,5-10  micron,  and  for  large  angle  scattering  (  ^  60  -  90®)  the 
effective  reflectivity  is  due  to  single  rather  than  to  multiple  scattering 
by  water  droplets,  as  a  result  of  the  relatively  large  absorption  and 
the  forward  peak  in  the  differential  scattering  cross  section.  The 
effective  diffuse  reflectivity  under  these  conditions  is  of  the  order  of 
1-2  percent.  For  actual  situations,  some  discussion  has  been  given  of 
atmospheric  transmission  losses  as  well  as  of  meteorological  data  on  the 
distribution  of  high  altitude  clouds.  Comparison  with  experiment  is  made 
in  a  (classified)  appendix,  and  the  results  are  found  to  be  consistent 
with  the  limited  experimental  data.  The  statistical  problem  consists  of 
the  detection  of  rare  high  radiance  signals  in  a  background  having 
occasional  high  radiance  levels.  The  detection  criterion  thus  involves 
maximizing  peak  rather  than  mean  signal-to-noise  ratios.  An  analysis 
in  terms  of  likelihood  ratios  is  developed,  and  it  is  shown  that  the 
Neyman-Pearson  criterion  provides  a  usable  partial  maximization  of  sig¬ 
nal  detection  probability.  In  practical  cases,  the  Neyman-Pearson 
criterion  is  equivalent  to  setting  up  a  single  threshold.  Limited  numeri¬ 
cal  applications  are  made  for  single  and  correlated  detector  elements. 
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SECTION  1 


INTRODUCTION 


1.1  THE  PROBLEM 

A  primary  factor  in  IR  surveillance  from  satellites  for  the 
purpose  of  detecting  rocket  boosters  is  the  infrared  background  in  which 
the  booster  must  be  detected.  This  report  deals  with  the  infrared  back¬ 
ground  in  the  2.7  micron  band  arising  from  clouds  illu^iinated  by  the  sun. 

It  is  known  that  sunlit  clouds  are  the  major  sources  of  background  radi¬ 
ance. 

The  problem  of  the  sunlit  cloud  background  can  be  separated  into 
two  parts.  The  first  part  is  the  determination  of  the  reflected  solar 
flux  from  an  illuminated  cloud  taking  properly  into  account  the  infrared 
absorption  in  the  atmosphere  itself.  The  second  part  is  the  specification 
of  the  statistical  character  of  the  background  in  terms  of  cloud  frequency 
distributions  for  height,  size,  and  for  irregularities.  Here  it  is  essential 
to  consider  various  kinds  of  detection  and  data  processing  systems  because 
different  statistical  properties  of  the  background  determine  the  perform¬ 
ance  of  the  different  systems. 

In  the  present  report  a  fairly  complete  treatment  of  the  first 
part  of  the  problem  is  presented.  The  second  part  is  also  treated  with  an 
attempt  to  achieve  some  generality.  However,  the  treatment  is  far  from 
complete  and  should  only  be  regarded  as  an  initial  attempt. 
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1.2  OUTLINE  OF  PRESENT  WORK 


Because  of  the  unsatisfactory  state  of  our  knowledge  of  the 
reflectivity  of  clouds  in  the  infrared  --  theoretical  estimates  differ 
by  almost  two  orders  of  magnitude,  a  thorough  scudy  of  infrared  reflec¬ 
tion  from  clouds  is  made  in  Sections  2  and  3,  The  problem  of  Mie  scat¬ 
tering  from  water  droplets  having  a  distribution  of  radii  is  treated  in 
Section  2,  where  it  is  shown  that  in  the  near  infrared  the  effective 
cloud  reflection  mechanism  for  large  scattering  angles  (0  ^  60°)  is 
single  scattering,  while  the  effective  cloud  reflectivity  for  these 
large  scattering  angles  is  of  order  1  -  2%.  Thus  the  problem  of  the 
scattering  of  infrared  from  clouds  of  water  droplets  has  been  here  solved 
with  sufficient  accuracy  for  most  practical  purposes,* 

The  results  of  Section  2  are  applied  in  Section  3  to  the  case 
of  a  real  atmosphere,  where  transmission  losses  and  scattering  effi¬ 
ciencies  are  combined  to  provide  a  series  of  isoradiance  plots  assuming 
uniform  tropospheric  cloud  cover.  Some  discussion  is  given  of  the  possible 
effects  of  large  ice  crystals  (Appendix  3A) ,  while  the  results  for  cloud 
reflectance  are  compared  with  the  limited  experimental  data  in  Appendix  3B, 
and  are  found  to  be  consistent  with  these  data. 


"k 

For  definiteness  the  numerical  analysis  is  limited  to  the  2.7  micron 
band.  The  same  scattering  analysis  with  somewhat  different  numerical 
coefficients  applies  for  a  2.7/^j^  A.  ^10  micron,  but  at  lower  wave¬ 
lengths  the  problem  is  essentially  different  because  the  absorptivity 
of  the  droplets  is  reduced  to  such  an  extent  that  multiple  scattering 
has  to  be  taken  into  account.  At  longer  wavelengths  it  is  essentially 
a  matter  of  grey  body  radiation. 
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The  treatment  of  the  statistical  aspects  of  the  background  in 
Section  4  begins  with  a  general  discussion  of  criteria  for  performance 
based  on  statistical  decision  theory.  Conceptual  difficulties  with 
elemental  costs  and  with  a  priori  probabilities  which  are  essential  con¬ 
cepts  in  this  formulation  are  pointed  out.  It  is  ultimately  made  plausible 
in  Section  4.2  that  a  simple  criterion  based  on  likelihood  ratios  has 
some  degree  of  validity.  The  likelihood  ratio  test  is  applied  to  various 
assumed  radiance  distributions  in  Section  4.4  where  it  is  shown  that  non- 
Gaussian  biraodal  distributions  have  an  entirely  different  effect  depending 
on  whether  the  system  is  evaluated  by  likelihood  ratio  or  r.m.s.  noise 
criteria.  In  such  bimodal  distributions,  the  large  background  radiances, 
although  infrequent,  have  a  large  effect  on  performance  which  is  not  re¬ 
flected  in  the  r.m.s.  noise  criterion.  It  is  also  shown  by  an  example 
in  Section  4.2  how  the  frequency  distribution  for  radiance  is  related  to 
the  height  distribution  of  the  clouds  through  the  calculation  of  cloud 
reflectivity  given  in  Sections  2  and  3.  In  Sections  4.3  and  4.4,  an 
attempt  is  made  to  apply  the  likelihood  ratio  test  to  a  comparison  of 
single  sensors  with  multiple  sensors  with  correlation,  leading  to  a 
rather  obvious  result  that  correlation  does  enhance  system  performance 
if  the  spatial  extent  of  a  cloud  is  larger  than  the  spatial  extent  of 
the  target.  Finally,  in  Section  4.6  some  remarks  are  made  about  optimum 
linear  filters. 

In  Appendix  4A  an  attempt  is  made  to  study  two-dimensional 
autocorrelation  functions  and  Wiener  spectra  on  the  surface  of  a  sphere. 

It  was  found  that  the  Wiener-Khinchine  theorem  does  not  hold  under  these 
conditions . 
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SECTION  2 


SCATTERING  OF  INFRARED  FROM  ASSEMBLIES  OF  WATER  DROPLETS 


2.1  PHYSICAL  PROPERTIES  OF  HIGH  ALTITUDE  CLOUDS 

A  high  altitude  cloud  is  an  assembly  of  water  droplets  or  ice 
crystals,  or  a  mixture  of  the  two  at  an  altitude  of  approximately  20,000 
feet  or  more.  Such  a  cloud  may  move  relatively  slowly  under  the  in¬ 
fluence  of  horizontal  winds,  as  in  the  case  of  cirrus  clouds  or  oro- 

3 

graphic  disturbances,  or  comparatively  rapidly  (10-10  cm/sec)  as  a 

2.1 

result  of  internal  convection  as  in  the  case  of  cumulus  clouds. 

As  yet  there  is  rather  little  detailed  experimental  data  on 
high  altitude  cloud  properties.  Measurements  of  low  altitude  clouds 
indicate  water  droplets  with  a  mean  radius  of  the  order  of  5  micron  with 

3 

perhaps  0.1-1  gm  H^O/m  .  At  higher  altitudes  one  may  perhaps  expect  a 
somewhat  lower  water  content.  The  droplet  size  distribution  may  be 
represented  as^*^ 


n(r)  dr  -  N^  r*^  e  dr  (2.1) 

where  n(r)  is  the  number  of  particles  per  cc  with  radius  r  -  (r  +  dr). 

2 . 2 

The  mode  radius  of  this  distribution  is  r  *  n/a.  Deirmendjian  *  uses 

-1  ^ 

n  “  6,  a  =  1.5  micron  ,  so  that  r  =  4  micron.  Here  this  distribution 

c 
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is  examined  in  detail.  Cases  with  smaller  and  larger  droplet  size, 

r  «  2,  8  micron,  and  with  the  same  mean  width  of  the  frequency  distri- 

^  -1 
butlon  as  for  n  =  6,  a  =  1.5  micron  ,  are  also  investigated,  as  are 

models  with  r  «  4  micron  but  with  narrower  and  wider  standard  deviations, 
c 

All  the  particle  size  distributions  are  shown  in  Figure  2.1  and  Table  2.1. 
In  general  we  shall  deal  with  thick  clouds--this  to  be  defined  precisely 
later“-and  thus  the  actual  number  of  particles  per  cc  or  the  liquid 
water  content  of  the  cloud  do  not  affect  the  scattering  directly. 

A  cloud  rises  due  to  the  buoyant  forces  present  when  the 
temperature  within  it  is  higher  than  the  ambient  temperature.  As  the 
cloud  rises,  its  temperature  drops  as  a  result  of  both  its  adiabatic  ex¬ 
pansion  and  of  turbulent  mixing  with  the  surrounding  atmosphere.  Thus, 
under  certain  conditions,  there  will  be  ice  crystals  in  clouds.  For 
instance,  cirrus  clouds  are  largely  or  entirely  ice,  while  cumulus 
clouds  may  also  contain  some  ice  crystals. 

There  are  wide  variations  in  the  shape  and  size  of  ice  crystals, 

ranging  from  frozen  water  droplets,  which  may  be  almost  spherical  and 

of  the  same  mass  as  the  water  drops,  all  the  way  to  prisms  and  hexagonal 

2. 1 

cylinders  up  to  500  micron  in  their  longest  dimension.  The  shape, 
size,  number  density  and  motion  of  these  ice  crystals  are  complicated  and 
incompletely  understood  functions  of  the  ambient  conditions,  and  it  is 
difficult  to  make  any  definitive  statement  about  them,  particularly  with 
respect  to  the  frequency  with  which  high  altitude  clouds  of  ice  crystals 
may  occur.  In  this  study,  clouds  will  be  treated  as  assemblies  of  water 
droplets.  In  other  words,  effects  arising  from  the  specular  reflection 


It  should  be  pointed  out  that  water  droplets  of  1-10  micron  radius 
freeze  at  -35  to  -40°C  because  of  surface  tension  effects.  See  Ref.  3.3. 
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TABLE  2.1 


SCATTERING  OF  INFRARED  FROM  VARIOUS  DISTRIBUTIONS  OF  PARTICLES 


Particle  Mode  Radius 

2 

4 

4 

4 

8 

»  n/a  (micron) 

Width  of  Distribution 

Regular 

Narrow 

Regular 

Wide 

Regular 

n 

2 

24 

6 

2 

24 

a  (micron 

1 

6 

1.5 

0.5 

3 

Mean  Mass  per  Particle 

2.5x10’^° 

3.40x10"^*^ 

-10 

6.25x10 

2.01xl0‘ 

2.72x10' 

(gni) 

3 

Density  of  Water  (gm/m 

)  .0251 

.0340 

.0625 

.2011 

.2723 

for  100  particles/cc 

O',  (10  ^  cm^) 

abs  '  ' 

.234 

.358 

.562 

1.247 

2.009 

O  ,,(10“^  cm^) 
scatt 

.735 

1.349 

1.466 

2.151 

3.064 

cr  /  cr' 

abs  scatt 

0.318 

0.265 

0.383 

0.580 

0.656 

c  c 

4.65 

9.3 

9.3 

9.3 

18.6 
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will  not 


of  light  from  ice  crystals  such  as  haloes,  sun  dogs,  etc,, 
be  taken  into  account.  Some  of  these  effects  are  discussed  in  Appendix  3A, 

The  refractive  index  of  water  in  the  infrared  has  been  tabu- 
2 . 5 

lated  by  Centeno,  *  The  results  in  the  2.7  micron  band  are  shown  in 
Figure  2,2,  There  is  no  information  on  the  temperature  variation  of  the 
refractive  index  of  water,  or  on  the  refractive  index  of  ice  in  this 
frequency  range.  The  most  noticeable  feature  is  the  fact  that  water 
absorbs  significantly  only  for  X  ^  2.7  micron.  The  physical  reason  for 
this  is  that  this  is  the  wavelength  corresponding  to  the  highest  molec¬ 
ular  vibration  frequency  of  the  H^O  molecule,  or  to  the  maximum  frequency 
of  the  spectrum  of  “lattice  vibrations”  of  water, 

2.2  SCATTERING  OF  INFRARED  BY  A  SINGLE  DROP 

The  scattering  of  electromagnetic  radiation  of  wavelength  X 
from  a  sphere  of  radius  r  is  characterized  by  the  parameter 

X  =  perimeter/wavelength  -  2  TC  r/  X  •  (2,2) 

If  X  1,  we  are  in  the  “geometrical  optics”  region,  while  Rayleigh 

scattering  holds  for  x  «  1,  However,  for  x  1,  which  is  the  case  in  the 

visible  and  near  Infrared  for  the  particle  sizes  under  consideration,  the 

7  6  7  7 

difficult  “Mie  Scattering”  analysis  applies.  *  *  *  To  solve  the  problem 
in  this  regime  it  is  necessary  to  solve  Maxwell's  equations  for  a  plane 
incident  wave  and  spherical  outgoing  scattered  wave  subject  to  the 
appropriate  boundary  conditions  on  the  scattering  sphere.  The  resulting 
expressions  giving  the  angular  distribution  and  total  scattering  cross 
section  are  given  in  Ref,  2,6,  p.  630  ff,  and  Ref,  2,7,  p.  114  ff.  In 
the  notation  of  the  latter  reference,  the  results  are 
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IK: 


Iliiilliiillillllillll 
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i 


(2.3a) 


(0)  =  I  (0)  I  ^  ;  i  (0)  =  I  S2  (0)  I  ^ 

00  „  ,  ,  r  a  p  ^(cos0)  1  /I 

,  .gv  _  7  2n  +  1  n  .11  - +  b  (d/d0)  P  \cos0) 

*1  =  1  n(n  +  1)  [  sin0  n  "  J 

n.  .  \  f  b  P  ^  (COS0)  1  I 

’  ('Q^  -7*^  ilnjtJi  _rL_S - +  a  (d/de)P  \cos0) 

^2  (9)  -  4  n(n  +  D  L  sin0  '  n  J 

n  “  1 


vj/„’(y)  ^n^y) 

'n  Vp^'(y)  ^^(x)  -  tn  Vf/^(y) 

m  s^^'(y)  V|/(x)  -  4"„(y) 
m  <4>'(y)  I^(x)  - 


1/2 


v|/^(z)  =  (7t  z/2)  +  1/2 


(z) 


X^(z>  ■  -(.Kzn)'-'^  i/jU) 


X„c»)  -  i-Kzm^'^  H„+  1/2'^’ <*> 


y  *  m  X 

m  ■  n'  -  in”  (note  the  minus  signi) 


(2.3b) 

(2.3c) 

(2.3d) 

(2.3e) 

(2.3f) 

(2.3g) 

(2.3h) 

(2.3i) 

(2.3J) 


2-7 


I 


T 

f 

2  1  I  ‘f  *  I  “«!  '} 

J-'  =  (2  7r  r^/x^)  V  (2n  +  1)  Re(a  +  b  )  (2.4b) 

ext  ^  n  =  1  '  “  “ 

where  i  (0)  and  ij^(0)  are  the  intensities  of  radiation  scattered  through 
angle  0  with  polarization  parallel  and  perpendicular  to  the  plane  of 

observation,  which  is  the  plane  defined  by  the  incident  and  scattered  rays. 

*  (2) 

P  (cos0)  are  associated  Legendre  polynomials,  J^(z),  (z)  and  H  '  ,  (z) 

n  y  y  V 

are  respectively  Bessel  functions  of  first  and  second  kind  and  Hankel 
functions  of  the  second  kind 


j,  (z)  =  J  ^  (z)  -  i  N  ^  (z) 


(2.5) 


m  is  the  complex  refractive  index  of  the  scatterer  defined  here  with  a 

negative  imaginary  part.  is  the  total  scattering  cross  section 

SC  o  c  c 

and 


cr  ^ 

ext 


(T'  + 

scatt 


(2.6) 


is  the  extinction  scattering  cross  section; 
section. 


is  the  absorption  cross 


These  expressions  have  to  be  evaluated  numerically,  using  a 
high  speed  computer.  For  the  present  purpose  the  computations  were 
carried  out  on  the  Aeronutronic  IBM  709  computer  using  a  program  locally 
called  CLOUD.  This  program  computes  the  differential  scattered  intensity 
per  particle  averaged  over  the  particle  size  distribution  (2.1)  in  units 
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2 

of  (2  7t'R/p^  )  ,  where  R  =  distance  from  scatterer  to  observer.  It 
applies  for  0.1.  The  intensity  shown  in  Figure  2.3  has  been  aver¬ 
aged  over  the  scattered  intensities  for  polarizations  parallel  and  perpen¬ 
dicular  to  the  plane  of  scattering 

i^^(9)  =  1/2  [  ill  (0)  +  (O)]  (2.7) 

The  program  automatically  computes  the  angular  distributions  i  ,  i,  , 

av 

i  and  the  total  and  scattering  cross  sections.  The  results  of  Figure  2.3 
and  Table  2.1  give  the  results  over  the  2,7  micron  band;  in  other  words, 
for  the  scattering  cross  section  of  Table  2.1,  for  example 

(2.725)  +  (2.775)] 

(2.8) 

where  the  complex  refractive  indices  shown  in  Figure  2.2  show  a  number 
of  interesting  characteristics. 

(1)  The  effect  of  averaging  over  a  particular  size  distribution 
is  to  smooth  out  the  rapid  variations  in  the  angular  scattering  distri¬ 
bution.  For  illustrations  see  Reference  2.2. 

(2)  The  angular  distribution  is  always  heavily  peaked  in  the 
forward  direction  and  the  ratio  of  forward  to  backward  scattering  is 
very  insensitive  to  the  details  of  the  distribution.  This  may  be  under 
stood  because  the  refractive  index  (n *  -  in**)  is  not  very  different  from 


ic  2  -1 

For  Deirmendjian' s  distribution  ""--that  is,  n  =  6,  a  =  1.5  micron 
in  Equation  (2.1)--90%  of  single  scattered  radiation  goes  into  angles 
of  less  than  36^,  while  the  scattering  through  angles  of  greater  than 
90^  is  essentially  isotropic. 
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Differential  scattering  cross  section 


Differential  scattering  cross  section 
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FIGURE  2.3  DIFFERENTIAL  SCATTERING  CROSS  SECTION  FOR  A  SPHERICAL  WATER  DROPLET 
AVERAGED  OVER  PARTICLE  SIZE  DISTRIBUTIONS  OF  EQUATION  (2.1). 


(b)  EFFECT  OF  VARHNG  THE  DISTRIBUTION  MODE  RADIUS  r  AT  CONSTANT 
DISTRIBUTION  WIDTH 


2-10  b 


one,  so  that  the  reflectivity  R  for  an  infinite  plane  slab 


2  2 

(n»  -  1)^  +  n»^ 

(n'  +  1)^  +  n"^ 


(2.9) 


is  quite  small,  actually  of  the  order  1  percent  (see  Figure  2,2).  In 
fact,  the  relative  amount  of  energy  scattered  through  angles  greater  than 
’TT/jl  -  J(  ;c/2)/J(0)  where 


r  IK 

J(0)  -  d«$  sin9  d0  i^^(0) 


(2.10) 


is  also  of  the  order  of  1  percent. 

(3)  The  different  particle  size  distributions  do  show  some 
difference  in  their  behavior  at  intermediate  angles  (10^:$  0  ^^120^), 
but  even  these  differences  are  not  very  large.  This  is  certainly  a 
consequence  of  the  smoothing  effect  of  the  particle  size  distribution. 

(4)  The  total  or  extinction  cross  section  of  Equation  (2.6) 

ext  2 

increases  with  Increasing  r^,  but  not  strictly  proportionally  to  r^  . 

This  effect  again  is  due  partly  to  the  effect  of  the  size  distribution, 
and  partly  to  the  oscillatory  behavior  of  P*  ^^9, 

Figure  13.14). 

(5)  It  is  noteworthy  that  <7^,  /  increases  with  in- 

aDs  scact 

creasing  r  >  even  though,  of  course,  no  change  in  refractive  Index 
c 

occurs.  This  may  be  understood  qualitatively  as  follows.  In  the  geo¬ 
metrical  optics  region  defined  by  x  ■  2  7tr/X^^  1>  for  a  sphere  of 

2 

radius  r  Is  proportional  to  the  surface  area;  that  is,  ^  oC  ^  • 

scare 

On  Che  other  hand,  Che  absorption  cross  section  Is  proportional  both  to 
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the  scattering  cross  section  and  to  the  path  length  within  the  absorber, 

3 

which  again  is  proportional  to  the  radius,  so  that  qC  r*  ,  It  is 

clear  from  this  why  one  might  expect  ^catt  increase  with 

increasing  r^  or  x^. 

(6)  The  small  maximum  in  i  (9)  at  0  =*  170^  is  a  residual 

2  2  2  4 

rainbow  *  ’  *  that  is  essentially  a  refraction  phenomenon  arising  from 

internal  reflection  combined  with  scattering.  Increasing  the  absorption 
will  reduce  the  importance  of  this  phenomenon, 

(7)  The  insensitivity  of  the  results  to  the  particle  size 
distribution  leads  one  to  suspect  a  similar  insensitivity  to  wavelength 
in  the  2  micron  -  15  micron  region,  since  the  refractive  index  changes 
comparatively  little. 

To  sum  up  the  results  of  the  Mie  scattering  calculation,  most 

of  the  scattering  from  a  single  drop  is  in  the  forward  direction,  and 

the  scattering  is  essentially  isotropic  for  scattering  angles  0  ^90^. 

The  result  that  the  effective  large  angle  reflectivity  is  of  the  order 

of  1  percent  seems  to  be  quite  unambiguous;  any  statement  to  the  contrary 

is  dubious.  Finally  /  -wO.3,  which  represents  a  very 

ao  s  sc  a  1 1 

significant  degree  of  absorption  for  the  multiple  scattering  problem. 

2,3  SCATTERING  BY  AN  ASSEMBLY  OF  DROPS 
a.  Introduction 

It  is  required  to  calculate  the  spectral  radiance  of  a  cloud  of 

water  droplets  under  conditions  of  given  (solar)  illumination.  That  is, 

2  ^ 

a  flux  of  radiation  F^  watt/cm  in  the  relevant  frequency  band  is  inci¬ 
dent  on  the  cloud,  and  we  ask  for  the  spectral  radiance  H  of  the  cloud. 


-4  2 

In  the  frequency  band  2.65  -  2.80  micron,  “  5.7  x  10  watt /cm  . 
(See  Ref.  2.8). 
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which  is  the  power  sent  into  unit  solid  angle  from  unit  surface  area. 

2 

Thus  H  is  measured  in  watt/cm  sterad.  Explicitly,  if  the  illuminated 
surface  area  of  the  cloud  is  the  radiance  may  be  expressed  in 

terms  of  the  differential  scattering  cross  section  per  unit  area,  Q^(i,0) 
of  the  cloud  as 


H(0,  4>)  Q^(i,  0)  (2,11) 

The  aim  here  is  to  evaluate  Q  (i,0). 

c 

In  this  section  it  is  first  shown  that  multiple  scattering  is 
unimportant  except  for  very  small  angle  scattering.  Then  a  general 
treatment  of  the  reflection  from  a  planar  cloud  is  given.  In  appendices 
we  give  a  model  to  make  quantitative  estimates  of  the  effect  of  multiple 
scattering,  and  also  a  discussion  of  the  scattering  from  spherical  clouds. 
The  discussion  of  spherical  clouds  has  been  relegated  to  an  appendix 
first,  because  it  is  not  as  important  as  the  case  of  planar  clouds,  and 
secondly,  because  it  has  not  proved  possible  to  give  a  numerical  solution 
of  the  same  degree  of  generality  as  for  planar  clouds. 

It  should  be  noted  that  the  results  of  the  present  section 
have  rather  wide  applicability  since  the  scattering  cross  section  i^^(0) 
is  essentially  independent  of  the  drop  size  distribution  over  a  fairly 
wide  range,  and  since,  in  addition,  the  result  for  a  large  cloud  is  in¬ 
dependent  of  the  number  density  of  droplets. 

b.  Single  vs.  Multiple  Scattering 

The  possible  effect  of  multiple  scattering  presents  considerable 
difficulties  in  the  case  of  a  large  cloud.  The  general  problem  of  radia¬ 
tive  transfer  with  partial  absorption  and  an  anisotropic  angular  dlstri- 

2.9  2  10 

bution  has  not  been  solved.  *  ^  ’  However,  in  the  present  situation 
of  predominantly  forward  scattering  with  non-zero  absorption,  the  effect 
of  scattering  of  orders  higher  than  single  scattering  has  been  shown  to 
be  small  except  for  very  small  scattering  angles.  (See  Ref,  2.9b). 
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This  result  can  be  understood  qualitatively  as  follows:  Most 

of  the  scattering  is  in  the  forward  direction.  For  Deirmendjian* s 
2  2 

distribution,  *  90  percent  of  the  single  scattered  energy  goes  into 

0  <  36°,  while  the  scattering  is  essentially  isotropic  for  9  >  90°. 

Also,  CT'  /  ^  0.3  (cf.  Table  2,1),  so  that  in  three  collisions 

^  abs  scatt 

the  total  intensity  is  down  to  1/e  times  the  original  intensity  due  to 
absorption.  By  this  time,  however  90  percent  of  the  remaining  energy 
would  be  scattered  into  angles  less  than'NiT'  •  36°  -  65°.  Thus, 

multiple  scattering  will  be  unimportant  for  scattering  through  angles 
greater  than  60  to  90°. 

2,9b 

This  same  result  has  been  obtained  by  Goldstein  by  an 

explicit  numerical  evaluation,  Goldstein  compares  first  and  second 

order  scattering  with  absorption  by  using  some  representative  Mie 

2.2 

scattering  results  of  Dermendjian.  *  In  this  way  he  finds  that  single 

scattering  far  outweighs  double  scattering  except  for  scattering  through 
★ 

very  small  angles. 

The  predominance  of  single  over  multiple  scattering  makes  the 
calculation  of  the  reflection  of  near  Infrared  radiation  from  clouds 
possible.  This  result  has  a  number  of  interesting  consequences: 

(1)  The  absorptivity  of  water  droplets  is  very  high  throughout 
the  region  ^  =  2.5  -  10  micron  and  above.  Thus,  the  present  analysis 
applies  throughout  this  frequency  range  with  appropriate  quantitative 
modifications. 

•k 

It  should  be  pointed  out  that  while  the  expansion  in  different  orders 
of  scattering  can  be  written  down  formally,  yet  for  anisotropic  angu¬ 
lar  distributions  at  least  the  real  importance  of  the  different 
orders  of  scattering  is  very  hard  to  estimate  except  by  explicit 
numerical  evaluation. 
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(2)  In  the  visible  region  the  absorptivity  of  water  is  very 
low.  Thus,  in  this  region,  a  large  cloud  may  be  represented  at  least 
roughly  as  a  diffuse  Lambert *s  Law  reflector,  because  of  the  large 
amount  of  multiple  scattering.  However,  the  justification  of  this  pro¬ 
cedure  is  by  no  means  trivial,  and  in  particular  we  do  not  know  how  to 
treat  the  region  1  micron  <x<  2.5  micron. 

(3)  In  the  far  IR  (^10  micron)  a  cloud  acts  as  a  grey  body 

emitter. 

(4)  It  follows  that,  except  possibly  in  the  visible,  the 
description  of  a  cloud  as  a  Lambertial  diffuse  reflector  is  a  very  poor 
one  on  physical  grounds.  Of  course,  it  may  turn  out  under  certain  con¬ 
ditions  that  the  concept  of  “diffuse  reflectivity**  may  be  useful  for  an 
empirical  description  of  cloud  reflection  if  the  effective  diffuse  re¬ 
flectivity  coefficient  does  not  vary  very  much  with  scattering  angle  in 
the  appropriate  region, 

(5)  For  low  altitude  clouds  there  will  be  a  significant  degree 
of  attenuation  due  to  the  absorption  of  water  vapor  in  the  cloud  in  the 
2.7  micron  band.  This  mechanism  is  not  effective  above  30  Kft,  where 
the  effective  extinction  of  the  beam  within  the  cloud  is  due  entirely  to 
absorption  by  the  droplets  rather  than  by  scattering.  Thus  if  there  are 
N  scatterers  per  unit  volume,  the  effective  extinction  distance  within 
the  cloud,  ^  is  given  by  the  relation 

(6)  The  extinction  distance  ^  ^  enables  a  precise  distinction 
to  be  made  between  small  and  large  clouds.  If  the  characteristic  distance 

scale  of  a  cloud  is  denoted  by  L,  then  for  L  A  one  is  dealing 

3  ^ 

with  a  small  cloud  composed  of  NL  droplets  all  scattering  independent¬ 
ly.  Thus,  the  total  scattering  cross  section  of  a  small  cloud  is  given  as 
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NL^  cr 


scatt 


(2.13) 


Q(small)^ 

scatt 


On  the  other  hand,  for  a  large  cloud  defined  in  terms  of  L  ^  |3  ,  the 

2  2 

effective  cloud  volume  seen  is  not  L  but  L  ,  since  the  incident 
beam  only  penetrates  ^  into  the  cloud.  Thus,  for  a  large  cloud 


(large) 

^(scatt) 


^  N 


,2 

JL  <y 

scatt 


scatt 
^  abs 


(2,14) 


In  other  words,  for  a  large  cloud  the  total  scattering  cross  section  is 

independent  of  the  number  of  scatterers  per  cc  or  of  the  scattering 

cross  section  per  scatterer  .  Of  course,  the  angular  distribution 

scatt 

does  depend  on 

— . 

abs 


c .  A  Planar  Cloud 

Consider  the  single  scattering  from  a  large  slab  of  material 
containing  N  scatterers  per  unit  volume  and  having  a  flat  top  surface 
with  an  illuminated  area  which  is  the  (xy)-plane  z  =  0  of  our 

system  of  coordinates  (see  Figure  2,4),  The  angles  of  incidence  and 
exit  are  respectively  i  and  e  measured  from  the  local  normal,  and  z  in¬ 
creases  going  into  the  scattering  medium, 

A  volume  of  element  dV^  has  a  total  scattering  cross  section 

N  dV.  ,  and  the  total  amount  of  single  scattered  radiation  sent 

A  scant 

into  unit  solid  angle  in  the  direction  (9,  i)  is 

=  F^N  ®  ^  (2.15a) 
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FIGURE  2,4.  GEOMETRY  FOR  SCATTERING  FROM  A  ] 
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[Sct/m],  -  Cr  i„(e)/J<0) 


(2.15b) 


where  J(0)  =  745  for  Deirmendjian* s  distribution  of  droplet  sizes  (n  =  6, 
a  =  1.5  micron  in  Equation  (2.1).  This  quantity  is  defined  in 
Equation  (2.10).  All  the  quantities  entering  here  have  been  defined 
either  in  the  text  of  in  Figure  2.4.  The  factor  e  gives 

the  effective  extinction  of  the  ray  which  is  scattered  through  an  angle  0 
at  the  point  A.  The  effective  extinction  is,  of  course,  due  to  absorption 
(see  Equation  (2.12)  and  the  discussion  immediately  preceding  it). 


irried  out  extending  J'J'  dx  dy 
over  and  letting  z  go  from  0  to  o4>  .  This  gives  the  result 


The  integration  over  dV  is  cai 

A 


(2.16a) 


f(i,e)  =  1/cos  i  +  1/cos  e 


(2.16b) 


Finally,  using  Equation  (2.12)  for^  ,  one  gets  the  result 


i  (e) 

o  ^scatt'  ^abs^  J(0)  f(i,e) 


(2.17) 


This  result  may  be  understood  physically.  For  a  large  cloud 

the  absolute  magnitude  of  the  single  scattered  energy  is  proportional  to 

the  incident  flux  F  and  to  the  illuminated  surface  area  S...,  but  does 
o  ill* 

not  depend  on  the  absolute  magnitude  of  the  scattering  cross  section  per 
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scatterer  or  on  the  number  of  scatterers  per  cc,  N.  The  angular 

scatt 

distribution  depends  on  the  differential  scattering  cross  section  from 
each  scatterer,  modified  by  the  extinction  within  the  cloud,  which  gives 
rise  to  an  “albedo'*  factor  (  /  )  and  to  a  geometrical  factor 

Lf (i,e)  ^  ,  See  item  (6)  of  Section  2,2. 


The  differential  scattering  cross  section  per  unit  illuminated 
surface  area,  Q^(i>  0)  of  Equation  (2,11),  is  plotted  in  Figure  2,5  for 
various  valuesof  i  in  the  important  case  of  coplanar  scattering  when  the 
incident  ray,  the  local  normal,  and  the  exit  ray  all  lie  in  the  same 
plane,  so  that 


i  +  e  +  0  =  7r 


(2.18a) 


f(i,  e)  =  1/cos  i  -l/cos(9  +  i)  (2.18b) 


It  is  possible  to  write  down  the  result  (2.17,  2.11)  in  terms 
of  Lambert's  Law  diffuse  scattering  with  an  effective  diffuse  reflection 
coefficient  ^  (i>  0)«  appropriate  expressions  are 


Q^(i,  0)  *  (i,  0)  /  IT  3  cos  i  cos  e 

(2.19a) 

C^o'/a/i^  0 

"*1  (i>  0)  ^  g  ^  f(i,  e) 

abs 

(2.19b) 
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Differential  scattering  cross  section  per  unit  illuminated 
.  V  /  O —  abs  \ 

area,  285  Q^(i,0).  (  285  =  { -  )  J(0)  ) 

^  I  ^  scatt  / 


FIGURE  2.5  SCATTERING  FROM  A  LARGE  PLANAR  CLOUD 


NOTE:  For  the  calculations  presented  here,  it  has  been  assumed 

that  the  incident  ray,  the  reflected  ray,  and  the  local 
normal  -  AI,  AE,  AK  of  Figure  2.4  -  all  lie  in  the  same 
plane. 
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Some  typical  numerical  values  of  7t  (i,  9)  are  given  here: 


i  =  e  =  0  ,  0  =  180°;  rj  =  .0099  (2.20a) 

i  =  e  =  5°  ,  0  =  170°;  =  .017  (2.20b) 

i  =  e  =  45°,  0  =  90°  ;  7^  *  .024  (2.20c) 

i  =  e  =  60°.  0  =  60°  ;  Tj  =  .060  (2.20d) 

i  =  0,  e  =  0  =  90°  >  ^  °  (2.20e) 


In  view  of  the  single  scattering  mechanism,  the  concept  of  a 
•'diffuse  reflectivity"  is  rather  unphysical.  Also,  the  angular  variation 
of  the  numerical  values  of  7^  (i,  0)  indicates  that  it  is  not  a  very  use¬ 
ful  concept  for  quantitative  discussion. 


These  values^  like  the  results  of  Figure  2.3»  refer  to  the  coplanar 
case,  so  that  they  give  maximum  values  for  the  scattering  cross 
section  or  reflectivity. 


2-21 


REFERENCE  LIST  FOR  SECTION  2 


2.1  B,  J.  Mason,  "The  Physics  of  Clouds",  Oxford  (1957), 

2.2  D,  Deirmendjian,  "Scattering  and  Polarization  of  Polydispersed 

Suspensions  with  Partial  Absorption",  Rand  Corporation,  RM-3228-PR 
(June  1962).  See  also  earlier  publications  by  the  same  author 
referenced  therein. 

2.3  F.  H,  Ludlam  6e  B.  J.  Mason, "The  Physics  of  Clouds",  Handbuch  der 
Physik,  Springer,  Berlin  (1957). 

2.4  M.  Minnaert,  "The  Nature  of  Light  and  Color  in  the  Open  Air", 

Dover,  New  York  (1954). 

2.5  M.  Centeno,  "Refractive  Index  of  Liquid  Water  in  the  Near  Infrared", 

J,  Opt.  Soc.  America  31,  244  (1941). 

2.6  M,  Born  6e  E,  Wolf,  "Principles  of  Optics",  Pergamon,  New  York  (1959). 

2.7  H.  C.  van  de  Hulst,  "Light  Scattering  by  Small  Particles",  Wiley, 

New  York  (1957). 

2.8  F.  S.  Johnson,  Editor,  "Satellite  Environment  Handbook",  p.  79, 
Stanford  (1961). 

2.9a  J.  S.  Goldstein,  "The  Infrared  Reflectivity  of  a  Planetary 
Atmosphere",  Astrophys.  J, ,  132,  473  (1960). 

2.9b  "Numerical  Solutions  of  the  Auxiliary  Equation  for  an  Inhomogeneous 
Planetary  Atmosphere",  talk  presented  at  the  11th  International 
Astrophysical  Symposium,  Liege,  July  1962. 

2.10  V.  R.  Stull  6e  G.  N.  Plass,  "Emisslvity  of  Dispersed  Carbon  Particles", 
J.  Opt.  Soc.  America^,  121  (1960). 


2-22 


APPENDIX  2A 


A  MODEL  TO  ESTIMATE  THE  EFFECT  OF  MULTIPLE  SCATTERING 


A  phenomenological  way  to  estimate  the  importance  of  multiple 
scattering  is  to  regard  the  cloud  as  made  up  of  a  single  scattering  layer 
one  extinction  depth  thick  and  a  diffuse  core  with  effective  diffuse  re¬ 
flectivity  7j  The  coefficient  *  is  obtained  by  matching  the  radiance 
given  by  this  model,  H*(0) ,  with  that  for  the  single  scattering  model  for 
a  given  value  9*  of  0. 

Instead  of  Equation  (2,16a)  we  have 


H*(9)  =  j  Q  dz 


=  H  (0)  f 


1  - 


cos  i  cos  e 
(2A.la) 

+  F^(  */7r  )  COS  i  cos  e 


where 


9* _ 

f(i,e)  cos  i  cos  e 

abs 


(2A.lb) 


i  +  e  +  0*  =  7r 

Clearly  the  choice  of  0*  is  not  unique.  0*  “TT  for  backward  scattering 
gives  i  =  e  =  0,  f(i,e)  =  2,  *  »  .0099  e”^  =  .00134.  The  function 

Q^*(i,0)  is  shown  in  Figure  2.5  for  i  “  0,  for  comparison  with  Q^. 
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It  should  be  noted  that  the  fact  that  Q  ^  Is  similar  to  Q 

c  c 

does  not  by  itself  prove  that  multiple  scattering  is  unimportant,  but 
merely  demonstrates  the  self-consistency  of  the  present  way  of  intro¬ 
ducing  multiple  scattering  in  terms  of  a  diffuse  core  underlying  a 
single  scattering  layer. 
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APPENDIX  2B 


SCATTERING  FROM  A  SPHERICAL  CLOUD 


It  has  not  been  possible  to  carry  through  the  analysis  of 
Section  2,  Paragraph  2,3c  for  scattering  from  a  planar  cloud  to  the 
case  of  a  spherical  cloud.  The  reason  for  this  is  that  we  have  not 
succeeded  in  evaluating  the  volume  integral  of  Equation  (2,15a)  for 
spherical  geometry.  In  the  present  appendix  are  given  the  formulation 
of  the  problem  and  an  approximate  estimate  of  the  cross  section. 

Consider  the  scattering  of  electromagnetic  radiation  through 
an  angle  0  by  a  sphere  of  radius  a.  In  particular,  a  typical  ray  is 
shown  in  Figure  2B.1.  The  ray  is  incident  on  the  sphere  whose  center 
is  at  the  origin  0  at  a  point  I  whose  spherical  polar  coordinates  are 
(a,  scattered  at  a  point  A(a-z,  and  leaves  the 

sphere  at  a  point  E(a,  6")*  The  angles  lOA,  AOE  are  denoted  by 

g  "  and  the  angles  0^,  d^A  are  denoted  by  i,  e  in  analogy  with 
* 

the  planar  case. 


It  should  be  noted  that  1  and  e  are  constant  for  the  planar  case  of 
scattering  through  a  constant  angle  0,  so  that  one  gets  different 
scattering  curves  for  different  values  of  i,  as  are  shown  in  Figure  2.5, 
This  is  not  the  case  for  spherical  geometry,  where  i,e  vary  for  all 
the  different  rays  corresponding  to  scattering  through  an  angle  0,  so 
that  there  is  a  unique  angular  scattering  curve.  However,  this 
feature  complicates  the  evaluation  of  the  integrals. 
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lAE  DENOTES  THE  PATH  OF  A  TYPICAL  RAY  SCATTERED  THROUGH  AN  ANGLE 
e :  NOTE  THAT  IN  GENERAL  THE  FOUR  POINTS  I,A,0,E,  DO  NOT  ALL  LIE  IN  A  PLANE. 
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The  triangle  rule  for  (lOA)  gives  the  exact  result 

^  “  ^rcsin  -  i  (2B.1) 

For  large  clouds  defined  in  analogy  in  the  discussion  of  Equation  (2,14) 
by 

,  ie  a»  z  ,  (2B.2) 

we  may  work  to  the  lowest  order  in  a/z.  Then  Equation  (2B.1)  gives 


t  '  *  (z/a)  tan  i  +  0(z/a)^  (2B.3a) 

r^  =  z/cos  i  +  0(z^/a)  ,  (2B.3b) 

and  to  this  order  the  volume  element  is 

=  dz  dS  [l  +  0(z/a)3  .  (2B.4) 

Here  dS  is  an  appropriately  defined  surface  element: 

dS  =  a^  g(  ,  5  )  d  co8^  d  $  (2B.5a) 

.  *  g'(’5'  8  ')  d  cosTy*'  d  8  '  ,  (2B.5b) 


where  g,  g*  may  be  functions  of  the  appropr5.ate  angles  but  are  Independent 
of  z  to  the  lowest  order  in  (z/a).  Thus,  to  this  order  the  radiance  of 
the  cloud  due  to  scattering  through  an  angle  9  is 
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(2b. 6a) 


[W/JA]  ,  -  N  [WwtJ,  J  IVj  »  ^ 


-  B(>’'/i'l3  ,  J  <IV, 


-  f(i,e) 


[l  +  0(z/a)j  (2B.6b) 


(2B.6c) 


where 


f(i,e)  =  1/cos  i  +  1/cos  e 


(2.16b) 


Thus  far  It  has  not  proved  possible  to  evaluate  the  Integral 
^  dS/f(l,e)  exactly.  In  general  one  may  put 


r _ ^ 

J  f(i. 


dS_  _  „  (0) 

e)  ■  ^ill 


<  l/f(i,e)^ 


av 


(2B.7a) 


2  f  />  ? 

Siii(0)  =  a  j  cos^/d  'I'  r  dU  =  20 


(2B.7b) 

where  the  coordinates  (  are  defined  in  Section  3. 3,  and  in 

particular  in  Figure  3.|  . 

An  upper  bound  to  ^l/f(i,e)^^^  may  be  set  as  follows.  First 
of  all,  consider  coplanar  scattering;  i.e.,  I,  A,  0,  and  E  all  in  the 
same  plane,  and 


l+e  +  0=7C-  e'-e"=7c£l-  0(z/a)][  (2B.8) 
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In  this  case  for  constant  scattering  angle  0, 


f(i,e)  »  1/cos  i  *  l/cos(9  +  i) 


(2B.9) 


and  as  i  varies. 


f(0)  *  2/sin(l/20)  ^  f(i,e)4fiP  (i  =  Tt /2)  . 


Thus  we  have  an  upper  bound  on 


<l/£(l.e)>^v 


<1/£(1,«))^^<  l/£(6)  . 


(2B.10) 


(2B.11) 


Thus  an  upper  bound  for  the  scattering  cross  section  per  unit 
illuminated  area,  given  by  the  expression 


%pher<») 


‘'o  ^ 


(Max) 


^spher 


(9) 


(2B.12) 


(Max) 

The  function  Qgpj^g^  (9)  is  shown  in  Figure  2.5. 

In  general  the  representation  of  a  planar  cloud  is  signifi¬ 
cantly  more  useful  than  that  of  a  spherical  cloud.  However,  under 
certain  conditions  of  oblique  viewing  the  effect  of  the  edges  of  clouds 
may  be  important,  and  here  the  model  of  a  spherical  cloud  may  prove  of 
interest.  In  any  case,  the  analysis  can  certainly  be  extended,  although 

numerical  techniques  may  be  necessary  to  evaluate  • 

^  ^  av 
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SECTION  3 


SCATTERING  OF  SOLAR  INFRARED  FRCW  CLOUDS 


3.1  METEOROLOGICAL  DATA 

Natural  clouds  are  made  up  of  assemblies  of  water  droplets 

or  ice  crystals.  In  either  case,  the  water  content  is  in  the  range 

3  3 

0.1-1  gm/m  at  the  higher  altitudes,  and  up  to  10  gm/m  at  sea  level. 

For  water  clouds,  observed  droplets  range  in  size  from  2-30  micron 

radius,  with  mode  radii  4-8  microns.  In  the  case  of  ice  crystals,  the 

particles  may  be  of  roughly  the  same  size,  or  there  may  be  significantly 

fewer  particles  up  to  500  microns  in  their  largest  dimension.  The 

shape  of  the  ice  crystals  may  be  hexagonal,  prismatic,  or  complex,  and 

the  crystals  may  possibly  be  aligned  by  convection  currents  within  a 

cumulus  cloud,  or  (less  probably)  by  wind  shear  in  the  case  of  cirrus  or 

stratus  clouds. 

Cumulus  clouds  in  general  cool  by  adiabatic  expansion  accompany¬ 
ing  their  rise,  rather  than  by  thermal  exchange  with  the  ambient  atmosphere. 
These  clouds  undergo  constant  convective  stirring,  which  leads  to  the 
characteristic  anvil  or  mushroom  shape  of  thunderheads.  It  should  be 
mentioned  that  1-10  micron  droplets  freeze  at  -35®C  to  -40^C  on  account 
of  surface  tension  effects,  so  that  the  question  of  under  what  conditions 
clouds  are  made  up  of  ice  crystals  rather  than  water  droplets  does  not 
necessarily  have  a  simple  answer. 


3-1 


Below  the  tropopause  (30-40  Kft)  there  is  a  great  deal  of 
local  atmospheric  structure  and,  hence,  various  types  of  clouds.  In 
general,  there  is  a  relatively  high  probability  (10%  for  the  continental 
U.S.A.)  of  more  or  less  uniform  cloud  cover  at  the  tropopause.  At  the 
present  there  is  not  much  quantitative  data  on  higher  clouds  except  for 
Tiros  satellite  pictures,  which  have  not  yet  been  processed  completely. 
Tropical  thunderstorms  associated  with  hurricanes  or  typhoons  form 
cumulus  systems  up  to  several  hundred  miles  across  which  have  been  known 
to  ’’burst  through*'  the  tropopause  and  rise  to  a  significant  height.  At 
higher  latitudes  there  is  not  much  statistical  information,  but  one  might 
infer  that  for  a  cloud  to  rise  a  height  h  above  the  tropopause  the  linear 
dimension  of  the  cloud  has  to  be  at  least  of  the  order  h. 

There  are  known  several  rare  and  not  very  dense  types  of  very 
high  altitude  clouds.  These  are  **noctilucent’' clouds  at  up  to  30  Km 
(90-100  Kft)  altitude,  which  are  very  thin  and  composed  of  ice  crystals, 
and  ’’nacreous"  or  "mother-of-pearl"  clouds  at  20-30  Km  which  are  inferred 
to  be  composed  of  spherical  particles  of  rather  uniform  size  on  account 
of  their  color. 

In  summary,  it  must  be  stressed  that  the  present  discussion 
clearly  does  not  do  justice  to  the  large  amount  of  meteorological  work 
of  observing  and  classifying  cloud  data.  Up  to  the  present,  most  experi¬ 
mental  work  deals  with  rather  low  altitude  clouds  (less  than  10  Km  or 
30  Kft)  as  far  as  observation  and  sampling  go.  New  information  from 
satellites,  balloons  and  high  altitude  jet  airplanes  is  beginning  to 
supplement  these  data.  Unquestionably,  in  3  to  5  years  there  will  be 
available  much  more  information  on  high  altitude  clouds  than  exists  at 
present. 
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3.2  TRANSMISSION  LOSSES 


a.  Intcoduction 

The  attenuation  of  solar  infrared  in  the  2.7  micron  atmospheric 
absorption  band,  which  is  here  interpreted  as  3675  +  100  cm“^,  is  to  be 
examined  for  heights  above  10  Km.  The  absorption  in  this  band  is  due  both 
to  CO2  ^nd  to  H2O.  The  CO2  in  the  upper  atmosphere  is  well  known  to  be 
mixed  uniformly  with  air  (.03%  by  volume).  The  amount  of  CO2  at  S.T.P. 
in  a  vertical  one-way  path  to  infinity  is  240  atm-cm  (cf.  Ref.  3.1, 
ch.  VIII),  and  at  a  pressure  of  (p/po)  atmospheres,  the  residual  amount 
is  240  (p/Pq)  atm- cm. 

The  amount  of  water  vapor  and  its  distribution  are  not  so 
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firmly  established.  Here  the  latest  compilation  of  Gutnick  *  is  used. 

The  results  are  shown  in  Table  3.1.  It  should  be  noted  that  Gutnick's 
estimates  only  extend  up  to  34  Km,  and  it  is  quite  possible  that  there 
is  a  significant  amount  of  water  at  higher  altitudes.  For  instance,  if 

o 

the  density  of  water  vapor  from  34-50  l&n  is  0.2  gm/m  ,  the  amounts  of 
precipitable  water  at  all  altitudes  in  Table  3.1  have  to  be  increased  by 
.0032  gm/cm  .  This  is  simply  a  measure  of  the  experimental  uncertainty 

'ic 

of  this  quantity. 


It  should  be  pointed  out  that  it  is  not  completely  clear  whether  the 

3.5 

water  vapor  concentration  is  a  single  function  of  height  *  ,  or  whether 

in  fact  there  are  large  variations,  by  factors  up  to  100  -  1000  between 

3  6 

a  **wet**  and  a  **dry”  atmosphere  *  .  The  present  concensus  favors  the 
uniform  distribution. 
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TABLE  3.1 


ATMOSPHERIC  WATER  VAPOR  CONCENTRATION  AT  HIGH  ALTITUDES 


Altitude  h 
(Km) 


Atmospheric  Density 
(Kg/m^) 


Mixing  Ratio 
m  (gm/Kg  air) 


Water  Density 
pm  (gm/m^) 


Precipitable 

Water 

(gm/cm^)  f£om 
34  Km  Down 


10 

4.135 

X 

o 

1 

.038 

1.571 

X 

10-2 

.00678 

12 

3.119 

X 

10-^ 

.017 

.530 

X 

10-2 

.00468 

14 

2.279 

X 

10-' 

.010 

.228 

X 

10-2 

.00392 

16 

1.665 

X 

IQ-' 

.0095 

.158 

X 

10-2 

.00353 

18 

1.217 

X 

10"' 

.012 

.146 

X 

10-2 

.00323 

20 

8.891 

X 

10-2 

.018 

.160 

X 

10-2 

.00292 

22 

6.500 

X 

10-2 

.027 

.176 

X 

10-2 

.00259 

24 

4.752 

X 

10-2 

.039 

.185 

X 

10-2 

.00223 

26 

3.436 

X 

10-2 

.060 

.206 

X 

10-2 

.00184 

28 

2.436 

X 

10-2 

.088 

.217 

X 

10-2 

.00141 

30 

1.786 

X 

10-2 

.125 

.223 

X 

10-2 

.00097 

32 

1.304 

X 

10-2 

.  180 

.235 

X 

10-2 

.00051 

34 

9.602 

X 

10-3 

.290 

.279 

X 

10-2 

0 

(1)  Precipitable  water  from  height  h^  to  h 


1  J  h 


(2)  If  ( ^  m)  s  2  X  lO’  gm/m  from  34  -  50  I&n  altitude,  the  additional 

2 

amount  of  precipitable  water  is  .0032  gm/cm  below  34  Ken. 
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The  transmissivity  for  a  vertical  one-way  path  from  infinity 
down  to  the  reference  altitude  h  is  given  in  Table  3.2.  The  transmissi¬ 
vity  factors  used  there  have  been  obtained  from  the  calculations  of 
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Plass,  Stull,  and  Wyatt  '  .  It  should  be  mentioned  that  the  transmissivity 
calculations  apply  to  constant  pressure;  we  have  used  values  corresponding 
to  \  p(h),  where  p(h)  is  the  atmospheric  pressure  corresponding  to  the 
reference  altitude.  This  is  clearly  not  completely  satisfactory  from 
a  quantitative  standpoint,  but  provides  an  estimate,  which  is  all  that  is 
called  for  in  the  present  application. 

b.  Slant  Path  Transmission 

The  fact  that  the  densities  of  CO2  and  H2O  are  functions  only 
of  altitude  (rather  than  of  geographic  location)  implies  that  for  a 
slant  path  from  a  point  A  to  infinity,  inclined  at  an  angle  to  the 
vertical,  the  effective  absorbing  path  length  is  given  by  the  relation 

Slant  Path  Length  (Angle  V  )  *  (Vertical  Path  Length  )  /  cos  (3.1) 

This  condition  holds  for  an  appropriate  range  of  2^  2  in  particular,  if 

=  radius  of  the  earth,  h  “  maximum  height  at  which  there  is  significant 
attenuation  (h^  ^  100  -  150  Kft) ,  then  (3.1)  holds,  provided 

(  I  -TT  )  »  \l  %  ^  (3-2) 

or,  since  ^  20  minutes  of  arc,  Eq.  (3.1)  holds  for  0  ^  V  <  39^. 

If  the  vertical  one-way  transmissivity  from  point  A  to  infinity 
is  t,  then,  provided  Eq.  (3.1)  holds, 

Slant  Path  Transmission  Loss  s  9  )  >  0  <  n  1  (3.3) 
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ty  factors  are  obtained  from  the  following  tabulated  values,  all 


u  3.8 
where 


n  -  1  for  the  "weak  line  approximation" 

n  =  ^  for  the  "strong  line  approximation  with  no  overlap"  (3.4) 
n  ^  ^  for  the  "strong  line  approximation  with  overlap" 
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The  overall  problem  is  a  very  complicated  one  *  .  In  the  calculations 
of  isoradiance  plots  in  Section  3.3,  we  use  n  =  1;  i.e.,  the  "weak  line 
approximation. " 

3.3  REFLECTION  OF  SOLAR  INFRARED  FROM  A  UNIFORM  CLOUD  COVER 
a.  Introduction 

In  Section  3.3  we  discuss  a  uniform  tropospheric  cloud  background. 
In  practice  there  would  be  large  black  areas  due  to  the  absence  of 
clouds  at  the  tropopause  and  occasional  bright  spots  arising  from  high 
altititude  clouds. 

For  definiteness  we  consider  a  satellite  D  at  a  height 
(r-a)  =  (  ^  “  l)a  above  the  surface  of  the  earth,  whose  radius  is  a. 

For  a  numerical  example,  we  take  the  satellite  height  to  be  2000  miles, 
so  that  ^  *=1.5.  Let  the  sun  S  be  at  an  angle  of  elevation  .  relative 

to  the  satellite.  The  question  is  now  the  following.  The  satellite 
scans  the  portion  of  the  surface  of  the  earth  that  it  can  observe  with  a 
very  narrow  beam.  Under  these  conditions,  how  much  radiation  in  the 
2.7  micron  band  does  the  satellite  receive  as  a  function  of  its  direction 
of  view,  of  the  elevation  angle,  of  the  normal  one-way  transmissivity 
factor  t  of  Section  3.2,  and  of  the  scattering  law  i^^(0)  of  Section  2 
and  Table  3.3? 
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TABLE  3.3 


THE  SCATTERING  FUNCTION  i  (9)* 

av 


Scattering  Angle  0  (degrees)  i^^(0) 

0  10,448 

5  5,970.2 

10  1,599.4 

15  559.9 

20  297.8 

30  107.6 

40  44.07 

50  20.10 

60  10.21 

70  5.929 

80  4.091 

90  2.659 

100  2.171 

120  2.078 

140  1.984 

160  1.891 

180  1.798 


This  function  differs  from  the  exact  result  of  Section  2  in  that  the 
“rainbow**  maximum  at  0  =  170^  has  been  smoothed  out. 
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The  numerical  calculations  have  been  carried  out  for  t  *  0.1  and 
0.5  and  for  oc  «  0^,  ^5°,  90^,  120°,  and  the  calculation  has  been  programmed 
so  that  a  different  scattering  law  inserted  readily  in  the 

form  of  a  table  if  it  is  desired  to  repeat  the  calculation  for  different 
values  of  the  parameters. 

b.  Geometrical  Considerations 

The  geocentric  coordinate  system  used  here  is  a  slight  modification 
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of  that  of  Zirker,  Whipple  and  Davis  .  All  the  coordinates  that  arise 
are  shown  in  Figure  3.2,  where  defines  the  z-axis  and  DOS  the  (zy) -plane, 
so  that 


D  «  (0,0, a  y  )  (3.5a) 

S  =  lim  t  (0,sin  oi  ,  cos  o<.  )  (3.5b) 

oO 

The  geocentric  angles  are  a  latitude  ^  and  a  longitude  oj  .  The 
latitude  Ip  =  0  defines  the  equatorial  (yz) -plane,  and  the  latitude 

defines  the  positive  x-axis.  The  longitude  CJ  -  0  defines  the 
positive  z-axis,  and  point  A  on  the  surface  of  the  earth  has  the  coordinates 

A  =  a(sin  cos  sin  GD  ,  cos  cos  cO  )  (3.6) 

and  the  surface  element  on  the  sphere  is 

dS  -  a^  cos  d  vj/  d(o  •  (3*7) 

The  vector  A?  has  magnitude  I),  and  is  given  as 

A^  e  a  (-sin  Vp  ,  -cos  \p  sin  6J  ,  ^  -cos  Vp  cos  )  (3.8a) 

/a^  ~  -  I  +  -  2  ^  cos  vp  cos  (O  (3.8b) 
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FIG.  3.2.  LIMITING  VALUES  OF  CO  . 

The  viewed  illuminated  portion  of  the  earth  is  characterized  by 

|u  (  ^  6  =  arc  )  because  of  the  condition  O^D  ^  iA„D  =  /r/2.  See 

Eq.  (3.12). 


In  terms  of  these  geocentric  angles  (  ,  CO  )  ,  the  local  angles 

of- incidence,  exit,  and  scattering  at  the  point  A  are  given  by  the  following 
relations: 


cos  i 

-  ^s> 

=  cos  ^ 

cos  (  o4.  -  (O  ) 

(3.9a) 

cos  e 

“  ^®0A’  ^AD^ 

~  (1/D)( 

^  cos  cos  OJ  -1) 

(3.9b) 

cos  ^ 

^  ^  cos  -  cos  ^  cos  (cK 

(3.9c) 

where  is  the  unit  vector  parallel  to  OA,  and  so  on.  The  angle  ^  =  TT  -  0, 
where  0  is  the  scattering  angle  of  Section  2, 

In  principle,  this  specifies  the  problem.  However,  it  is  more 
convenient  to  express  the  results  in  terms  of  the  angles  of  view  of  the 
satellite,  here  chosen  as  (  ^  )  and  sketched  in  Fig.  3.1.  By  referring 

to  Fig.  3.1,  one  sees  that 


cos  COA  =  ~  ^cos  CO 

(3.10a) 

CA  =  a  sin  c3a  =  a(l-cos^V|/  cos^OJ 

(3.10b) 

AB  =  a  sin  , 

(3.10c) 

sin  ^  »  AB/AC  =  sin  Vp  /  (l-cos^\p  cos^CO 

(3.11a) 

sin  >1  *  AC/AD  *  (l-cos^Vp  cos ^  CO  )^/j4. 

(3.11b) 
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FIG.  3.1.  GEOMETRY  FOR  SCATTERING  OF  SOLAR  RADIATION  FROM  A  POINT  A  ON  THE  SURFACE 
OF  THE  EARTH  TO  THE  DETECTOR  D. 

cX  is  the  elevation  angle,  (U  ,ti‘)  are  the  latitude  and  longitude  of  A,  i,  e,  0  = 

'ft  -  are  the  angles  of  incidence,  exit  and  scattering  for  the  ray  SAD  under 
consideration  here.  ,  J  )  the  angles  of  view  of  A  from  the  satellite. 


There  is  one  further  matter  to  be  discussed,  namely  the  limiting 
values  of  the  angles  (  vJ/,CO).  First  consider  60  ,  and  see  Fig.  3.2. 

The  limiting  ray  DA^  corresponds  to  DA^O  =  ‘TC/a.  ,  and  thus  for  0  ^ 

Min  (S‘,JoC-  ^  S  =  arc  cos  (1/  ^  ) .  (3.12) 

For  the  limit  on  ,  see  Fig.  3.3,  where  the  limiting  ray  DA^  is  shown. 

cos  DOA^  =  1/^  =  (eJq,  =  cos  y/  cos  (JJ  (3 . 13) 


so  that  for  each  CO  which  satisfies  (3.12),  we  have  the  limits  on 


^  arc  cos(l/  0  cos  CO  ) 


(3.14) 


c.  Isoradiance  Plots 

The  earth  is  taken  to  be  covered  by  a  uniform  layer  of  thick 
flat-topped  clouds  at  an  altitude  corresponding  to  a  vertical  one-way 
transmissivity,  t,  measured  from  the  cloud  top  altitude  to  infinity.  One 
calculates  the  spectral  radiance  H(A)  at  a  variable  point  A  on  the  layer 


H(A) 


r  T  i 

^  scatt  av 

L  O-abs  J 


(0) 


.f(i.e) 


J(0)  f(i,e) 


(3.15) 


The  point  A  is  characterized  by  geocentric  latitude  and  longitude  ^ , CO 
respectively  and  for  a  given  set  of  angles  (  vf,  OJ  )  the  angles  of  Incidence 
exit  and  scattering  i,e,  0  •  TT -  are  given  by  Eq.  (3.9).  The  term 
t^(i»®)  gives  the  fractional  energy  loss  on  a  slant  path:  sun-cloud  point  A 
detector  under  the  "weak  line  approximation"  of  Eq.  (3.3,  3.4).  In  general, 
this  factor  is 


f  (i.e) 


.  ^l^cos  i)"  +  l/(cos  e)^  ^  0  ^  ^  ^ 


(3.16) 
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FIG.  3 . 3.  LIMITING  VALUES  OF  LATITUDE  '  •  . 

For  a  given  value  of  ,  the  viewed  illuminated  portion  of  the  earth  is 
characterized  by  O^D  >■  O^dD  =  /2,  and  thus|Y|  ^  arc  cos(l/  /'  cos  See 

Eq.  (3.12-14). 
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2  3.10 

e  .000568  watt/cm  is  the  solar  radiance  in  the  2.7  micron  band  ’  , 

while  the  remaining  terms  just  give  the  angular  scattering  cross  section 

per  unit  area,  Q^(i,0)  of  Eq.  (2.11,  16,  17)  and  Fig.  2.5. 

The  scattering  function  i^^(0)  used  here  is  given  in  Table  3.3. 
It  differs  slightly  from  the  accurate  form  of  Section  2.2,  Fig.  2.3,  in 
that  the  *‘rainbow'*  maximum  at  0  =  170^  has  been  smoothed  out.  This  was 
done  because  this  maximum  complicates  the  plotting  significantly,  and  its 
effect  is  a  factor  of  two  at  most,  which  is  not  particularly  significant 
for  a  first  look  at  the  isoradiance  contours. 

The  computational  problem  is  the  following.  For  a  given  cloud 
point  characterized  by  the  geocentric  latitude  and  longitude  ( Cp  ,  ) > 

the  satellite  viewing  angles  (  ^  determined  from  Eq.  (3.11)  and 

the  angles  (i,e,  ^  =7r  '  obtained  from  Eq.  (3.9).  The  function 

i^^(0)  is  obtained  by  a  *'table  look-up*'  routine  and  H(A)  calculated  from 
Eq.  (3.15).  This  part  of  the  problem  is  carried  out  on  a  high-speed 
digital  computer  for  a  range  of  points  A.  From  these  data  isoradiance 
contours  are  obtained  as  a  function  of  the  satellite  look  angles  (  ^  ) 

by  manual  cross  plotting.  The  isoradiance  plots  for  a  range  of  trans¬ 
missivity  factors  t  and  elevation  angles  oC,  are  given  in  Fig.  3.4. 


3-15 


Fig.  3.4  Isoradiance  Plots 

The  lines  drawn  are  contours  of  equal  brightness  as  a  function 
of  the  satellite  view  angles  (  ^  ^  They  are  calculated  from  Eq.  (3.15) 

for  various  elevation  angles^  and  transmission  factors  t.  The  numbers 
listed  give  log  where  H  is  measured  in  watt/cm  steradi  increasing  up 
to  the  value  of  the  central  maximum  on  a  given  plot.  Graphs  are  shown  for 
-  0,  45°,  90°,  120°  and  for  t  -  0.1,  05. 
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FIGURE  3.4  (Continued)  ISORADIANCE  PLOTS 
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FIGURE  3.4  (Continued)  ISORDIANCE  PLOTS 
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FIGURE  3.4  (Continued)  ISORDIANCE  PLOTS 
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FIGURE  3.4  (Continued)  ISORDIANCE  PLOTS 
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FIGURE  3.4  (Continued)  ISORDIANCE  PLOTS 
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FIGURE  3  4  (Continued)  ISORDIANCE  PLOTS 
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Effect  of  Inhomogeneities 

Apart  from  the  obvious  effect  of  holes  in  the  cloud  cover  at 
the  tropopause  on  the  isoradiance  curves  of  Fig.  3.4,  irregularities  and 
detailed  structure  in  the  top  of  the  cloud  cover  can  produce  certain 
changes.  Casual  observation  from  jet  airplanes  shows  that  sometimes  the 
cloud  cover  looks  extremely  uniform,  sometimes  it  has  a  very  regular 
structure,  and  sometimes  it  is  quite  irregular,  with  convective  humps 
above  the  general  level  of  the  surface  and  with  occasional  wisps  of  ice 
crystals  from  cirrus.  Of  course,  by  looking  from  airplanes  one  cannot 
give  a  quantitative  measure  of  the  scale  of  irregularities.  The  question 
arises,  how  much  and  in  what  ways  will  these  phenomena  affect  the  general 
level  of  the  radiance  of  the  cloud  cover? 

First  of  all,  the  relevant  measure  of  ^'flatness"  of  the  cloud 
tops  is  given  in  terms  of  the  extinction  length 


For  Deirmendjian's  distribution  of  Section  2  (cf.  Ref.  2.2,  Table  2.1) 
with  N  =  100  particles/cc ,  the  extinction  distance  is  of  the  order  of 
200  meters. 

Consider  next  regular  modulations  of  the  cloud  top.  If  the 
cloud  top  is  made  us  of  closely  packed  circular  cylinders  of  spheres, 
elementary  geometrical  considerations  show  that  if  the  sun  local  normal 
angle  is  less  than  60°,  the  effect  of  shadowing  will  be  to  cut  the  over¬ 
all  radiance  by  less  than  a  factor  of  two.  The  effect  of  modulation 
clearly  depends  on  the  frequency  response  of  the  detector  system  as  well 
as  on  the  precise  direction  of  the  scan,  but  no  detailed  examination 
seems  called  for  at  this  point. 
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The  effect  of  large  irregularities  will  clearly  be  to  distort 
the  picture  significantly,  but  any  estimate  of  the  magnitude  and  character¬ 
istics  of  these  distortions  involves  so  much  detailed  input  information 
that  once  again  no  further  investigation  seems  appropriate  in  the  present 
context . 

Finally,  it  is  clear  that  the  overall  effect  of  inhomogeneities 
depends  on  the  resolving  power  of  the  detecting  system. 
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APPENDIX  3A 


SPECULAR  REFLECTION  FROM  ICE  CRYSTALS 


At  present  there  exists  essentially  no  experimental  information 
on  the  actual  composition  of  high  altitude  clouds  ,  but  it  seems  at  least 
plausible  to  suppose  that  they  contain  a  certain  fraction  of  ice  crystals. 

It  may  well  be  that  the  reflection  properties  of  natural  clouds  of  ice 
crystals  do  not  differ  very  much  from  those  of  clouds  of  water  droplets, 
but  this  point  is  certainly  not  established.  The  purpose  of  this  appendix 
is  to  outline  some  characteristic  differences  between  water  droplets  and 
ice  crystals,  mainly  to  provide  a  starting  point  for  more  detailed  future 
work. 

The  problems  to  be  mentioned  are  the  composition  as  far  as 
crystal  size  under  given  conditions  is  concerned,  the  motion  of  the  crystals 
with  reference  to  possible  alignment  effects,  and  the  diffraction  pattern 
produced  by  crystals  of  given  size. 

b.  Composition  and  Structure  of  Clouds 

It  is  established  that  the  size  and  crystal  structure  of  ice 
crystals  depends  significantly  on  the  conditions  under  which  a  cloud  is 
produced. 

(1)  Current  work  by  B.  J.  Mason*s  group  at  Imperial  College, 

3,11 

London,  has  established  both  experimentally  and  theoretically  that 
cloud  droplets  of  less  than  20  fX  radius  cannot  capture  smaller  droplets 
to  grow  by  coalescence.  On  the  other  hand,  in  a  supercooled  mixture  of 
water  droplets  and  ice  crystals,  the  crystals  grow  at  the  expense  of  the 
droplets  as  a  result  of  a  vapor  pressure  gradient. 
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(2)  Artificial  ice  fogs  produced  in  winter  by  seeding  the 

3.12 

moist  air  near  Old  Faithful  geyser  in  Yellowstone  National  Park  appear 

to  show  an  overall  angular  scattering  distribution  rather  similar  to  that 

3. 13 

predicted  by  Mie  theory  calculations  like  Deirmendjian's  and  the  present 
ones  (cf.  Section  2).  However,  there  are  also  some  very  striking  specular 
effects . 


(3)  As  against  this,  certain  other  measurements  on  natural 

ice  fogs  off  the  West  Coast  of  Alaska  dealt  essentially  with  frozen  water 

j  1  3.14 

droplets 

3.2 

(4)  Extensive  measurements  by  Weickmann  6e  Aufm  Kampe  have 
discovered  prismatic  single  crystals  up  to  several  hundred  microns  in 
length  and  with  a  length-breadth  ratio  of  the  order  of  1  -  5,  and  hexagonal 
plates  of  the  same  maximum  dimension  and  perhaps  50  microns  thick. 

In  summary,  there  does  not  appear  at  present  to  be  sufficient 
experimental  evidence  to  enable  us  to  make  any  categorical  statements  about 
the  size  and  shape  of  ice  crystals.  In  particular,  there  seems  to  be  no 
justification  in  favor  of  using  results  for  low  altitude  ice  fogs  to  infer 
the  behavior  of  high  altitude  cirrus  clouds.  Such  factors  as  the  absolute 
humidity,  density,  temperature  and  temperature  gradient  and  availability 
and  type  of  condensation  nuclei  are  simply  not  understood  adequately  as 
far  as  their  detailed  effect  on  cloud  composition  and  structure  is  concerned 


b.  Motion  and  Alignment  of  Ice  Crystals 

The  terminal  fall  speed  of  a  sphere  of  radius  a  and  density  ^  in 


a  medium  of  density 


m 


and  viscosity 


is 


m 


term 


=  (2/9)  8  a  (  ^  -  ^  ^)/  >1 


(A.l) 
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For  air,  -yi  =10  gm/cra  sec,  so  that  for  water  droplets  of  radius  a  *  1, 

»  ^  -2  2 
10,  100  microns  the  terminal  speed  is  respectively  10  ,  1 ,  10  cm/sec. 

In  other  words,  for  crystals  of  dimensions  1  -  500  microns,  possible 

3 

convective  flow  speeds  in  clouds,  which  may  go  up  to  10  cm/sec,  are  adequate 
to  maintain  the  particles  in  suspension. 

Ice  crystals  may  be  aligned  either  by  non-uniform  flow  fields 
such  as  wind  shear  (especially  for  needle-shaped  particles),  or  by  aero¬ 
dynamic  forces,  which  are  particularly  important  for  disks  and  other 
flattened  objects  (L/D  >  1).  The  effect  of  wind  shear  can  readily  be 

shown  to  be  unimportant.  In  laminar  flow,  wind  shears  ^  u/  ^  y  are  of 

3  15 

the  order  .03  (cm/sec)/cm  *  ,  and  thus  for  a  crystal  of  dimensions 

10^  X  10^  X  10  (in  micron),  the  force  A  rl  u/ y  is  of  order  10"*^^  <iyn, 

^  —  4 

as  against  the  gravitational  force  mg  ~  10  dyn.  Thus,  unless  velocity 
gradients  in  turbulent  flows  are  10^  -  10^  times  larger  than  these  observed 
values,  the  effect  of  wind  shears  will  be  negligible. 

As  far  as  the  overal  aerodynamic  motion  of  objects  of  various 
shapes  under  the  action  of  viscosity  and  gravity  is  concerned,  qualitatively 
this  is  a  very  complex  *'phugoid‘'  motion.  Two  examples  of  this  are  shown 
in  Fig.  3A.1.  Whether  the  motion  is  damped  or  autorotating  depends  on  para¬ 
meters  such  as  the  wing  loading,  Reynold *s  number,  ratio  of  the  moments 

of  inertia,  etc.  The  overall  aerodynamic  problem  is  a  very  difficult 
3.16,  17,  18 

one 

c.  Diffraction  Pattern 

As  a  measure  of  the  importance  of  diffraction  effects,  consider 
the  Fraunhofer  diffraction  from  a  circular  aperture  of  radius  a.  The 
angle  0^  of  the  first  zero  in  the  diffraction  pattern  produced  by  light 
of  wavelength  X  is  given  by  the  relation 

sin  0^  =  0.61  (X  /a).  (A. 2) 
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(a)  DAMPED  MOTION  (b)  AUTOROTATING  MOTION 
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FIGURE  3A.1.  POSSIBLE  PHUGOID  MOTIONS  OF  ICE  CRYSTALS 


Clearly  this  angle  0^  gives  a  measure  of  the  effective  width  of  the 
diffraction  pattern.  For  /a)  ^  1,  0^  is  so  large  that  the  structure 
of  the  scattering  object  is  unimportant:  in  this  situation  treating  the 
crystal  as  a  sphere  will  clearly  be  a  good  approximation.  On  the  other 
hand,  diffraction  effects  are  clearly  unimportant  for  (  X  /a)  1,  so 

that  the  specular  reflection  will  be  important  in  this  region.  In  this 
geometrical  optics  case  it  is  necessary  to  understand  the  orientation 
and  motion  of  the  crystals,  as  well  as  their  shape  and  size,  in  consider¬ 
able  detail.  This  sort  of  information  is  simply  not  available  at  present. 

d.  The  Energy  Reflected  Specularly  from  a  Cloud  of  Ice  Crystals 

Let  the  flux  of  solar  radiation  be  F  watt/cm^  in  the  relevant 

o 

frequency  range.  Under  conditions  of  specular  reflection,  a  fraction 


q  = 


scat  t 


nr  + 
^scatt 


cr 


abs 


(A.3) 


of  this  energy  goes  into  a  solid  angle  ScJ  about  the  direction  of 
specular  reflection,  where 


fco  =  +  ScO 

diffraction  motion 


Sco 

diffraction 


ScO 


motion 


~  <  >  »v 


(A.  4a) 
(A.  4b) 

(A. 4c) 


where  @2  is  the  mean  change  in  orientation  of  the  crystal  under  the  action 
of  aerodynamic  forces. 
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For  a  single  crystal  of  area  A  =  TT  a  ,  the  energy  radiated 
per  unit  solid  angle  within  a  cone  about  the  direction  of  specular 

reflection  is 


-F^Aq/Sfio  .  (A.  5) 

Now  if  there  are  N  crystals  per  cc,  the  extinction  distance  is  of  the 
order 


(A.  6) 


Thus  for  a  large  cloud  the  total  energy  radiated  specularly  per  unit  of 
illuminated  surface  area  and  per  unit  solid  angle  within  a  mean  solid 
angle  about  the  direction  of  specular  reflection  is 


H 

spec 


_ 


11 


"ill 


(A..  7) 

2 

(watt/cm  sterad), 


where  q 


0. 7  while  Sui  Is  of  course  not  known  in  general, 


The  result  (A. 7)  should  be  compared  with  the  result  (2.17)  of 
the  Hie  scattering  analysis 


H  (0) 


F 

o 


cr 

scatt 

^bs 


i 

av 


(0) 


J(0)  f(i,e) 


(2.17) 


For  scattering  angles  greater  than  60°  or  so,  i  ^  (0)/f(i,e)  -w  1, 

av 

and  thus 


H  (0)  ^  F^/300  (watt/cm^  sterad)  for  0  ^60°-90°.  (A. 8) 


3-31 


It  is  clear  that  at  least  for  scattering  angles  greater  than  60^  -  90°, 
at  the  appropriate  angle  the  effect  of  specular  reflection  will  predominate 
over  Mie  scattering.  The  relative  magnitude  and  thus  the  lower  limit  to 
scattering  angles  for  which  this  is  true  depends  on  StO  $  but  for 

1,  0  X  60  90°,  H  ^  200  H  (9). 

’  spec  ' 
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SECTION  4 


DETECTION  AND  DECISION  SYSTEMS  FOR  NON-GAUSSIAN  BACKGROUNDS 


4 . 1  INTRODUCTION 

A  central  question  encountered  in  the  design  of  IR  surveillance 
detection  and  decision  systems  is  the  discrimination  between  a  target  and 
a  severe  background.  There  are  really  two  related  problems,  namely: 

(1)  How  to  find  the  optimum  detection  system. 

(2)  How  to  best  describe  the  background  and  target  radiances. 

For  our  purposes  discrimination  between  a  target  and  a  severe 
background  means  that  the  detection  system  will  be  required  to  decide 
whether  a  particular  alarm  signal  received  is  due  to  a  target  or  a  back¬ 
ground  condition,  that  is,  the  system  must  be  able  to  distinguish  between 
targets  and  severe  backgrounds  occurring  separately  rather  than  simul¬ 
taneously.  The  occurrence  of  a  target  together  with  a  severe  background 
will  be  a  very  rare  event  (and  therefore  relatively  unimportant)  since 
each  individual  event  in  itself  is  rare.  It  is  the  events  which  produce 
false  target  indications  which  cause  the  main  difficulty  in  IR  surveillance. 

For  example,  consider  an  infrared  detection  system  operating 
from  an  Earth  satellite  whose  purpose  is  to  detect  a  missile  launch  from 
plume  radiation.  If  it  is  accepted  that  there  exist  occasional  severe 
backgrounds  due  to  high-altitude,  infrared-reflecting  clouds,  then  the 


y 
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next  question  is  how  to  make  a  detection  system  work  despite  the  severe 
background.  Numerous  discrimination  schemes  have  been  suggested  which 
take  advantage  of  background  and  target  signal  characteristic  differences. 
Each  scheme  may  be  examined  in  terms  of  the  background  model (s)  which 
may  be  developed  to  decide  which  appear  most  promising.  Statistical 
methods  are  generally  involved,  and  such  methods  will  be  used  for  selec¬ 
tion  of  background  descriptions  and  comparison  of  various  alternate 
detection  systems. 

An  optimum  detection  system  is  defined  for  our  purposes  as 
one  which  is  capable  of  extracting  a  target  signal  from  a  background  and 
of  deciding  for  or  against  alarm  activation  with  a  minimum  of  incurred 
errors  or  '*cost."  There  are  two  kinds  of  possible  errors  which  may  be 
made  in  deciding  for  or  against  an  alarm  with  associated  probabilities 
of  error.  These  are: 

(1)  Mistaking  a  target  for  a  cloud. 

(2)  Mistaking  a  cloud  for  a  target. 

The  second  type  of  error  becomes  serious  if  the  cloud  radiance  probability 
density  distribution  overlaps  the  target  radiance  distribution  and  if 
the  threshold  for  alarm  decision  is  determined  so  that  a  warning  will  be 
sounded  for  the  majority  of  possible  target  radiances.  False  alarm  then 
becomes  a  nuisance,  and  it  is  imperative  that  the  detection  system  be 
designed  with  an  optimum  capability  for  discrimination  against  background 
radiances  in  the  high  radiance  tail  of  the  probability  density  distribution. 

A  general  criterion  for  system  performance  based  upon  statisti¬ 
cal  decision  theory  is  formulated  in  Section  4.2.  This  provides  a  funda¬ 
mental  basis  for  examining  various  discrimination  schemes.  The  general 
approach  has  been  to  minimize  the  expected  costs  for  the  various  possible 
detection  outcomes.  This  approach,  however,  requires  a  knowledge  of  the 
a  priori  probability  of  target  presence  and  the  assignment  of  costs  to 
the  possible  detection  outcomes. 

'fc 

See  Section  4.2  b  for  specification  of  the  four  possible  detection 
outcomes . 
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It  is  shown  in  Section  4.2  that  the  Neyman-Pearson  criterion^ 
of  minimizing  the  conditional  probability  of  target  miss  for  a  given 
conditional  probability  of  false  alarm  represents  a  partial  minimization 
of  the  expected  cost.  The  virtue  of  this  criterion  is  that  it  does  not 
require  a  knowledge  of  a  priori  probabilities  and  costs.  Furthermore, 
in  most  of  the  situations  encountered  in  the  present  study,  this  criterion 
reduces  to  the  single  signal  threshold  criterion  for  deciding  for  or 
against  target  presence. 

In  Section  4.3  the  properties  of  probability  density  functions 
representing  the  background  radiance  patterns  under  the  various  imposed 
conditions  (sun  angle,  elevation  angle  of  azimuthal  scan,  etc.)  are 
investigated  and  typical  non-Gaussian  models  are  presented.  Finally, 
since  the  physical  model  developed  permits  only  the  derivation  of  proba¬ 
bility  density  functions  for  background  radiance  and  not  for  the  signal 
input  to  the  decision  device,  the  conditions  under  which  these  two  proba¬ 
bility  density  functions  are  identical  are  examined,  and  possible  approxi¬ 
mate  methods  for  transforming  the  radiance  probability  density  function 
to  signal  density  function  are  considered. 

It  may  not  be  possible  to  discover  the  optimum  detection  system 
design.  In  fact,  one  would  probably  be  content  with  a  system  which  works 
despite  severe  background  without  too  much  regard  for  whether  the  system 
is  optimum  or  not.  A  variety  of  design  approaches  are  possible  if  optimi¬ 
zation  is  not  required  but,  instead,  simply  a  workable  system; 

A  performance  criterion  developed  in  Section  4.2  e  is  applied 
in  Sections  4.4  and  4.5  to  the  cases  of  the  single  sensor  and  the  multiple 
sensor  with  correlation  of  output  signals.  The  criterion  permits  the 
comparison  of  the  relative  performance  of  these  two  types  of  systems.  The 
result  of  the  comparison  indicates  that  correlation  can  enhance  system 
performance  considerably  if  sunlit  clouds  are  indeed  several  times  larger 
than  the  target. 
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In  Section  4*6,  the  optimum  linear  filter  is  considered.  In 
this  connection,  the  best  representation  of  background  infoarmation  has 
become  controversial.  As  Robinson^  has  pointed  out,  many  of  the  possible 
representations  of  background  radiance  information  do  not  retain  enough 
of  the  background  information  to  permit  the  design  of  an  optimum  detection 
system.  The  problem  of  representing  the  background  adequately  becomes 
more  severe  when  sophisticated  discrimination  methods  are  considered. 

Some  of  the  possible  representations  of  background  which  may  often  be 
found  inadequate  are: 

(1)  One  dimensional  Wiener  transform  and  autocorrelation  function. 

(2)  Two  dimensional  Wiener  transform  and  autocorrelation  function.^ 

(3)  Radiance  probability  density  function. 

(4)  Radiance  joint  probability  density  function.^ 

(5)  Rough  radiance  maps  using  limited  gray  scale. 

The  problem  with  the  use  of  the  first  two  descriptions  which 
Robinson^  has  pointed  out,  arises  because  of  two  factors.  First,  the 
background  probability  density  function  is  generally  non-Gaussian.  In 
this  case  the  methods  applied  to  the  selection  of  an  optimum  linear  signal 
processing  do  not  permit  the  selection  of  the  best  possible  processing 
which  may  actually  be  a  non-linear  filter  rather  than  linear.  Secondly, 
the  requirement  of  decision  for  or  against  alarm  of  the  basis  of  signal 
level  relative  to  a  threshold  introduces  an  inherent  non-linearity  in  the 
system.  Any  design  approach  which  does  not  account  for  the  non-Gaussian 
nature,  of  the  background  and  the  non-linearity  of  the  decision  processes 
cannot  be  expected  to  succeed. 

It  is  clear  that  it  is  not  essential  to  retain  all  background 
information  for  design  purposes,  since,  depending  upon  the  detection  and 
discrimination  method  selected,  much  information  will  be  irrelevant  to 
system  performance.  A  complete  background  description  would  consist  of 
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an  enormous  quantity  of  high  resolution,  large  radiance-ranges,  radiance 
maps  with  geographical  location,  time  (including  season  and  sun  angle), 
weather  conditions,  detector  position,  wavelength  region,  and  other  condi¬ 
tions  assigned  to  each  map.  Fortunately,  the  complete  description  will 
probably  never  be  required  for  the  design  of  any  particular  system.  The 
problem  is,  then,  reduced  to  one  of  selecting  the  appropriate  background 
description  for  the  design  of  a  detection  system  employing  the  selected 
discrimination  method. 

The  comparison  criterion  developed  and  applied  in  the  later 
sections  specifies  the  best  way  to  express  the  background  in  the  appropri¬ 
ate  form  with  regard  to  the  particular  discrimination  method  which  it  is 
desired  to  use.  Thus,  the  problem  of  the  appropriate  form  of  the  back¬ 
ground  representation  disappears.  However,  the  great  practical  difficulty 
of  obtaining  the  numerical  values  for  parameters  appearing  in  the  background 
expression  still  remains.  The  solution  to  this  particular  problem  lies 
in  improving  methods  of  data  gathering,  processing,  and  reduction  and  in 
performing  additional  background  measurements  of  sufficient  extent. 

4.2  PERFORMANCE  CRITERIA 

Various  performance  criteria  will  be  discussed  and  compared  in 
this  section.  The  general  approach  to  the  design  of  optimum  detection 
systems  has  been  to  formulate  the  expected  cost  in  terms  of  the  decision 
outcomes.  This  approach  is  one  suggested  by  statistical  decision  theory.^* ^ 
The  criterion  for  this  general  approach  is  that  the  best  system  is  one 
which  minimizes  the  expected  cost  of  the  system  outcomes.  Here  cost  is 
used  in  a  very  broad  sense  of  damage  or  hurt  of  any  sort  incurred.  How¬ 
ever,  the  practical  feasibility  of  this  decision  theory  approach  is 
limited  by  the  knowledge  of  certain  doubtfully  obtainable  a  priori  proba¬ 
bilities  and  by  the  necessity  of  assigning  cost  values  to  particular  detec¬ 
tion  system  outcomes.  We  will  see  that  a  partial  minimization  of  costs 
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is  possible  even  though  the  numerical  values  of  the  a  priori  probabilities 
and  costs  of  system  outcomes  are  not  known.  This  partial  minimization 
is  adopted  in  the  subsequent  applications. 

The  minimization  problem  is  outlined  as  follows:  A  sample 
function  is  observed  over  a  finite  time  interval.  This  sample  function 
consists  of  background  signal  having  known  statistics  and  may  or  may  not 
contain  a  target  signal.  The  observed  sample  function  is  suitably  pro- 
cesb  d  and  then  a  decision  is  made  as  to  whether  a  target  signal  is  present 
or  not.  For  each  of  the  four  possible  combinations*  of  target .signal 
presence  and  decision  a  cost  is  assigned.  The  system  is  then  optimized 
by  finding  the  decision  rule  which  leads  to  a  minimum  value  of  the  expected 
cost. 

a.  Probability  Density  Functions 

It  is  convenient  to  formulate  the  detection  problem  on  a  discrete 
basis.  The  background  is  described  by  a  set  of  numbers  which  give  the 
radiance  in  each  resolution  cell  of  a  finite  area.  Since  these  numbers 
are  obtained  by  sampling  the  outputs  of  sensors,  they  also  include  all 
of  the  noise  components  arising  from  the  background  signal  and  the  system. 

A  convenient  method  of  representing  these  background  samples  is  a  row 
matrix  of  n  elements,  where  n  is  the  number  of  resolution  cells  in  the 
area  being  considered.  The  components  of  this  matrix  may  be  assigned  in 
any  known  manner  to  the  resolution  cells.  It  is  assumed  that  the  time 
sequence  in  which  the  cells  are  observed  precludes  the  possibility  of  a 
particular  moving  target  appearing  in  more  than  one  cell.  Let  the  back¬ 
ground  matrix  be  designated  as 

'  ■  {'>l-  >>2 . 


* 


See  Section  4,2  b 
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The  background  has  associated  with  it  a  probability  density  function 
designated  as  w^(B)  and  normalized  such  that 

jw^(B)  di  =  1  (4.2) 

The  function  w^(B)  represents  the  probability  that  the  background  signal 
will  be  found  between  B  and  B  +  dB  when  it  is  multiplied  by  dB. 

The  target,  or  targets,  may  also  be  described  by  a  matrix, 
most  of  whose  elements  are  zero.  It  is  assumed  that  each  target  is 
sufficiently  small,  so  that  it  can  occupy  only  one  resolution  cell.  This 
matrix  is  designated  as 

^  '  {'^l’  ^2’ . ’  *^n} 

It  is  also  necessary  to  assign  a  probability  of  target  presence 
and  a  density  function  for  target  amplitude  when  a  target  is  present. 
Hence,  let  be  the  a  priori  probability  that  no  target  is  present,  and 
be  the  a  priori  probability  that  one  or  more  targets  are  present.  As 
will  be  seen  later,  it  is  often  difficult  to  say  precisely  what  these 
probabilities  are.  Since  these  probabilities  are  mutually  exclusive 

P,  +  P  =1  (4.4) 

1  o 

The  probability  density  function  for  target  amplitude  is  designated  as 
Wj^(T)  and  normalized  so  that 

jw^(T)  dT  =  1  (4.5) 

Thus,  the  a  priori  probability  density  function  associated  with  the 
target  matrix  is 

p(T)  =  S(T-O)  +  (T)  (4.6) 
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where 

/  -- 

1  for  T  »  0 

6(T-0)  =  j  (4.7) 

I  0  for  T  =  0 

The  observed  resultant  signal  is  in  general  some  combination 

of  the  background  signal  and  the  target  signal.  The  background  and 

'fc 

target  signals  are  assumed  to  be  statistically  independent.  For  most 
systtuis  the  target  and  background  signals  are  combined  by  simple  addition, 
but  it  is  not  necessary  to  make  this  assumption  for  the  general  derivation 
being  considered. 

Let  the  observed  resultant  signal  be  designated  by  a  row  matrix 
also,  so  that 

i  =  ^s^,  S2.  .  .  ..  S^J  (4.8) 


There  are  conceivable  conditions  where  background  and  target  signals 
are  not  statistically  independent.  For  example,  launching  might  be 
difficult  under  weather  conditions  that  produce  high  backgrounds. 

‘k 

If  the  target  is  a  true  point  source,  the  assumption  of  additivity  is 
exact.  Any  real  target  will,  however,  obscure  a  portion  of  the  back¬ 
ground.  Thus,  a  substitutive  rule  will  be  used  in  a  later  section 
since  the  target  image  for  systems  of  interest  will  be  on  the  order  of 
a  resolution  element  in  size.  This  approximation  also  permits  a  simpli¬ 
fication  in  the  mathematics. 
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Associated  with  this  matrix  are  two  conditional  probability  density 
functions,  p  (S  ) T  ^  0)  and  p  (S  I  T  =  0) .  The  first  probability  density 
function  is  the  probability  that  if  a  target  signal  is  present  and  has 
a  particular  value  T  then  the  background  signal  will  be  such  as  to  give 
a  resultant  signal  which  lies  between  S  and  S  +  dS.  The  second  proba¬ 
bility  density  function,  namely  p  (S  |  T  =  0) ,  is  the  probability  that 
if  a  target  signal  is  absent  then  the  resultant  signal  lies  between 
S  and  S  +  d'S.  This  latter  density  function  is  identical  with  the  proba¬ 
bility  density  function  of  the  background,  w^('S). 

Although  the  probability  density  functions  introduced  in  the 
preceding  paragraph  appear  simple  at  first  sight,  they  contain  a  wealth 
of  information  from  the  2n  random  variables  of  background  and  target* 

The  practical  description  of  the  background  in  terms  of  such  probabilities 
may  not  be  possible  at  all  as  will  be  seen  in  Section  4.3  where  this 
subject  is  again  taken  up.  Despite  the  complexity  of  these  conditional 
probability  density  functions,  they  enable  the  definition  of  various 
kinds  of  averages  which  are  practical  to  use.  Also,  restriction  of  the 
number  of  variables  permits  the  practical  applications  of  these  functions. 
In  any  case,  this  framework,  while  complicated,  is  necessary  to  describe 
the  problem  carefully  and  does  yield  suitable  quantities  for  practical  use. 


b ,  Decision  Rules  and  Cost  Functions 

The  decision  operation,  which  decides  for  or  against  alarm 
activation,  is  designated  by  ^  0  ) .  The  matrix  S  is  the  measure¬ 

ment  space  and  is  some  selected  subset.  For  a  single  threshold 
decision,  the  specification  of  Q  divides  the  measurement  space  into  two 
regions.  All  S  not  contained  in  Q  lead  to  a  decision  of  no  target 
present,  and  S  lying  within  Q  lead  to  a  decision  that  a  target  signal 
is  present.  Thus 

[  1  for  S  c  0 

(4.9) 

0  for  S  0  Q  =0 


S(S,  Q  )  = 
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The  sub-space  over  which  we  have  elected  to  decide  for  target 
presence  may  be  found  by  integration  of  S (^>  0  )  over  all  S.  Similarly, 
the  sub-space  over  which  we  decide  for  target  absence  is  the  integral  of 
{  1  -  5(^.Q  over  all  S. 

Finally,  it  is  necessary  to  assign  evaluation  or  cost  and  gain 
functions  to  the  four  possible  outcomes  which  are  the  result  of  time, 
target  condition  and  the  applied  decision  rule.  The  possible  outcomes 
and  heir  associated  cost  functions  C  (T,  S)  are: 

Decision  for  target  absent  when  target  is  absent.  ^ 

C  (T  ==  0,  5=  0)  =  C 

o 

Decision  for  target  present  when  target  is  actually  absent.  I 

C  (T  =  0,  5=  1)  -  I 

Decision  for  target  absent  when  target  is  actually  present  (4.10) 

c  (f  0,  5=0)  = 

Decision  for  target  present  when  it  is  present.  / 

C  (T  ^  0,  6=  1)  =  / 

The  costs  associated  with  the  correct  decisions  are  C  and 

o 

where  is  the  cost  of  deciding  target  absence  correctly  and 
is  the  cost  of  deciding  for  target  presence  correctly.  It  will  be  seen 
that  these  particular  costs  are  actually  gains  and  so  will  be  negative 
numbers  or  zero.  The  costs  associated  with  incorrect  decisions  are 
positive  costs  and  are  C£  and  C^,  costs  of  false  alarm  and  of  target 
miss,  respectively.  The  costs  C£  and  will  usually  be  assigned  large 
values  relative  to  the  magnitudes  of  Cq. 


4-10 


Since  the  occurrence  of  various  cost  values  is  random,  it  is 


to  find  the 

expe 

cted  value 

of  the  cost,  E  (T,  5)^ 

• 

<c  (T, 

6  )> 

=  E 

<Co>  + 

E  <^0 

(4.11) 

ly  seen 

that 

r 

❖o> 

=  C 

0 

1 

J 

1  P  ( s  1  T  = 

0) 

fl  -  S(S.Q  )}  dS 

(4.12) 

<C£> 

p 

P  (S  1  T 

=  0) 

5(S,  Q)  ds 

(4.13) 

=  c 

n 

1  p  (s  |t 

^  0) 

(l  '  6(S,Q)}  ds 

(4.14) 

<Cd> 

=  Cd 

f 

P  (six 

0) 

6(S,  Q  )  dS 

(4.15) 

from  which 


C=PiC  +P  C  + 
1  m  0  0 


«5(S,  n  )  |t  =  0)  (C^-C^) 

+  p  (S  I  T  0)  (C^  -  C^)j 


(4.16) 


dS 


In  deriving  Eq.  (4.16)  the  normalizing  relationships 

/ 

)  (S  |t  =  0)  dS  =  1 


and 


p  (S  j  T  =  0)  d^:  =  1 


(4.17) 


(4.18) 


have  been  used.  The  expected  cost  may  also  be  written  in  terms  of  the 
conditional  probability  of  false  alarm,  P£.  and  the  conditional  probability 
of  detection,  P^,  defined  respectively  by 


6(s,n  )  p  (S  (x  =0) 


dS 


(4.19) 
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and 


S(s,  0  )  P  (S  |t  ^0)  ds 


In  terms  of  and  the  expected  cost  is 


C=P,  C  +P  C  +P  P-  (C.-C  )  +  P,  P ,  (C,-C  ) 
Im  oo  of  fo  Id  dm 


(4.20) 


(4.21) 


The  general  design  criterion  may  now  be  stated.  The  optimum 
system  will  be  that  system  whose  parameters  produce  a  minimum  for  the 
expected  cost.  In  other  words,  the  decision  rule  S(s,n  )  must  be 
chosen  so  as  to  minimize  the  expected  cost.  The  appropriate  decision 
rule  is  that  S  (S,  Q)  zero  when  the  integrand  of  Eq.  (4.16)  is  positive 
and  unity  when  the  integrand  is  negative,  that  is 


Sis,  Q  )  = 


0  when 


1  when 


y\> 


C  -  c 

f  O 

c  -  c , 

m  d 


c.  -  c 

f  o 

c  -  c , 

m  d 


(4.22) 


where  the  likelihood  ratio  ,  is  defined  as 


yV  =  2l3..  0)  (4.23) 

P(S  I  T  =  0) 

The  actual  application  of  Eqs.  (4.22)  and  (4.23)  clearly  depends  upon 
having  some  data  with  respect  to  cost  values  and  a  priori  probabilities. 

Thus,  this  result  has  somewhat  limited  practical  scope.  One  way  to  get 
around  this  is  to  use  a  partial  minimization  of  the  expected  cost 
Eq.  (4.21).  This  leads  naturally  to  the  Ne3niian-Pearson  criterion  if  we 
minimize  the  last  term  on  the  right  of  Eq.  (4.21)  as  we  shall  see. 
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According  to  the  Neyman-Pearson  criterion,  a  decision  for 
target  presence  is  made  when  the  ratio  yv  defined  by  Eq.  (4.23)  exceeds 
some  particular  value,  say  K,  for  the  observed  S.  When  for  some  S,  y\. 
falls  below  K,  then  a  decision  for  target  absence  is  made.  The  value  of 
K  is  selected  to  give  a  preassigned  false  alarm  rate.  This  amounts  to 
a  specification  of  Cf.  For  the  Neyman-Pearson  decision  rule,  the  decision 
operator  is  then 

^  yv  >K 


(s,  n )  =  < 


1  when 


0  when 


(4.24) 


K 


A  unique  decision  operator  is  not  specified  since  K  may  have  any  desired 
value. 


It  is  of  interest  to  compare  Eqs.  (4.22)  and  (4.23)  with  the 
Neyman-Pearson  criterion  for  selecting  the  decision  rule.  In  the  previous 
case  a  definite  value  of  A  is  specified  when  the  costs  and  a  priori 
probabilities  are  known.  Under  the  Neyman-Pearson  criterion,  however, 
the  threshold  is  arbitrary. 

The  following  argument  indicates  that  the  Neyman-Pearson 
criterion  is  a  first  step  towards  minimizing  the  expected  cost.  If  the 
expected  cost  is  expressed  as  in  Eq.  (4.21),  a  possible  approach  to  mini¬ 
mizing  the  cost  would  be  to  optimize  one  of  the  four  terms  on  the  right 
side  while  holding  the  other  three  constant.  The  only  control  permitted 
is  over  the  terms  involving  Pf  and  P<j.  Under  the  Neyman-Pearson  criterion, 
the  choice  of  K  fixes  the  value  of  Pf.  Since  always  a  negative 

quantity,  a  partially  optimum  system  will  be  that  for  which  P^  is  a  maxi¬ 
mum  so  that  the  cost  is  minimized.  This  procedure  corresponds  precisely 
to  the  Neyman-Pearson  criterion. 
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The  great  virtue  of  the  Neyman-Pearson  criterion  lies  in  the 
fact  that  it  enables  some  control  over  the  cost  without  requiring  a  speci¬ 
fication  of  the  costs  or  a  priori  probabilities.  This  advantage  is  to 
some  extent  illusory,  however,  since  it  is  not  possible  using  this  criterion 
to  know  what  the  actual  cost  is.  Furthermore,  since  costs  are  not  evaluated, 
the  actual  gain  achieved  by  using  the  optimum  system  rather  than  some 
compromise  system  cannot  be  determined. 


c .  Single  Threshold  Criterion 

Another  decision  criterion  often  used  is  the  single  signal  threshold 
rule.  Under  this  rule  a  decision  for  target  presence  is  made  when  ‘S’  exceeds 
some  threshold  value,  say  S^.  The  decision  operator  in  this  case  is 

f  1  when  S  ^  "Sq 


Q)  -  { 


(4.25) 


0  when  S  ^ 

This  criterion  is  the  same  as  the  Neyman-Pearson  criterion  when  yv  is 
a  monotonic  function  of  S  in  which  case  the  Neyman-Pearson  condition 
separates  the  signal  space  into  just  two  regions.  However,  if  for  a  given 
value  of  y\  there  correspond  several  values  of  S,  then  no  single  signal 
threshold  can  satisfy  the  Neyman-Pearson  minimization. 


An  example  will  be  discussed  to  demonstrate  the  relation  between 
the  Neyman-Pearson  condition  and  the  single  threshold  condition.  Suppose 
that  the  background  probability  density  function,  Wq(5)  is  bimodal  with 
most  clouds  providing  radiance  less  than  the  target  radiance  and  a  very 
few  clouds  producing  radiances  much  in  excess  of  the  target  radiance  values. 
Note  that  these  latter  radiances  correspond  to  small  values  of  the  proba¬ 
bility  density  function.  The  Neyman-Pearson  condition  will  require  two 
signal  thresholds:  one  which  separates  the  low  cloud  radiances  from  the 
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range  of  target  radiances  and  a  second  one  which  separates  the  high  cloud 
radiances  from  the  target  range.  In  this  manner  both  low  and  excessively 
high  cloud  signals  will  be  excluded  from  the  ^  subspace.  However,  the 
region  of  high  radiances  which  is  excluded  by  the  upper  signal  threshold 
will  not  affect  system  performance  significantly  provided  the  probability 
of  occurrence  of  such  clouds  is  very  small,  i.e.,  their  contribution  to 
p  (S  It  =  0)  d"?  is  small.  For  all  practical  purposes,  then,  a  single 
signal  threshold  would  suffice.  This  is  expected  to  be  the  case  for 
missile-launch  detection  systems.  Thus,  the  establishment  of  a  single 
signal  threshold  for  practical  systems  is  equivalent  to  design  according 
to  the  Neyman-Pearson  criterion  provided  that  the  target  and  cloud  distri¬ 
butions  are  well  separated  and  most  backgrounds  are  in  one  region  of  the 
signal  space. 

Other  decision  criteria  are  possible  such  as  for  systems  using 
correlation.  For  example,  a  decision  for  target  presence  might  be  made 
only  when  the  threshold  for  a  single  detector  is  exceeded  and  the  signals 
from  adjacent  detectors  remain  below  the  signal  threshold.  In  this  way 
discrimination  against  large  reflecting  clouds  can  be  achieved. 

If  the  detection  system  is  equipped  with  data  storage  capability, 
the  possibility  of  reserving  a  decision  until  more  data  is  accumulated  is 
permitted. 

d.  A  Priori  Probabilities  and  Costs 

As  already  mentioned,  the  Neyman-Pearson  approach  is  the  only 
approach  which  is  reasonably  general  yet  independent  of  a  priori  proba¬ 
bilities  and  cost  assignments.  In  order  to  use  the  more  general  minimum 
cost  criterion  of  Eq.  (4,22)  these  probabilities  and  costs  must  be  known 
or  estimated.  It  might  prove  instructive  to  examine  what  is  involved  in 
estimating  these  quantities. 


In  this  case,  Wald*s  sequential  criterion  might  be  useful,^ 
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Since  the  probability  of  any  enemy  ICBM  attack  varies  from  day 
to  day  as  the  international  situation  changes,  the  estimates  of  must 
be  updated  since  is  the  a  priori  probability  of  target  presence  in  the 
interval  t^  and  t^  +  dt.  Other  conditions  which  are  usually  steady 
relative  to  the  height  of  internal  tension  feed  into  the  estimates  of  P^^ 
such  as  the  number  of  enemy  missiles  on  launch  pads  and  their  probable 
targets.  Various  governmental,  military,  and  intelligence  agencies  make 
estimates  of  the  probability  of  an  enemy  attack  from  which  it  might  be 
possible  to  derive  the  probability  of  target  presence  in  the  field  of 
view  of  the  detection  system,  Pj^.  Since  much  subjectivity  goes  into  these 
estimates,  their  reliability  may  be  poor.  So  far  no  test  cases  have 
occurred  to  verify  the  estimates,  although  past  experience  and  the  absence 
of  an  attack  indicate  that  P^^  is  usually  considerably  smaller  than  unity. 
In  view  of  the  questionable  reliability  of  estimates  of  P^^,  any  design 
criterion  which  is  independent  of  a  knowledge  of  P^  is  favored. 

The  cost  functions  of  Eq.  (4.10)  are  even  more  difficult  to 
estimate.  In  any  discussion  of  costs  of  errors  of  a  surveillance  system, 
complicated  interactions  of  the  surveillance  system  with  other  parts  of 
the  defense  system  have  to  be  considered  which  causes  difficulty  in  such 
a  general  discussion  unless  the  rules  of  the  defense  game  are  known.  What 
is  the  cost  of  a  detection  miss  for  one  single  element  in  an  extensive 
surveillance  system  consisting  of  warning  radars  as  well  as  satellite 
infrared  systems?  Perhaps  the  cost  is  that  of  not  having  the  advantage 
of  an  earliest  possible  warning.  Thus,  it  might  be  inversely  proportional 
to  the  extent  of  build-up  of  the  retaliatory  forces  actually  achieved  at 
the  moment  the  attacking  warheads  reach  their  targets.  Or  perhaps  it  is 
a  function  of  the  warning  time  actually  achieved  by  the  overall  system 
so  that  the  cost  of  a  miss  depends  upon  the  functioning  of  the  rest  of 
the  surveillance  system. 
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What  is  the  cost  of  a  false  alarm?  Does  it  amount  to  the  cost 
of  scrambling  SAC  and  initiating  the  procedure  for  additional  retaliation? 
Does  it  also  depend  upon  the  lack  of  credibility  in  the  warning  provided 
by  the  rest  of  the  surveillance  system? 

One  might  be  permitted  to  say  that 

if  the  assumption  that  retaliation  is  not  actually  ordered  until  one  or 
more  enemy  hits  have  been  experienced  is  allowed.  In  this  case,  a  false 
alarm  will  initiate  an  alert  and  also  will  be  checked  out  by  future  experi¬ 
ence,  but  a  miss  may  mean  the  loss  of  the  advantage  of  an  alert  and  possibly 
of  some  of  the  retaliatory  forces  as  well. 

So  far  the  decision  criteria  for  an  optimum  system  and  for  other 
systems  have  been  discussed.  The  assumption  which  is  made  is  that  the 
system  which  uses  the  minimum  cost  decision  rule  would  be  the  optimum  system. 
Since  it  is  difficult,  if  not  impossible,  to  know  and  the  costs,  it 
will  not  always  be  possible  to  discover  the  optimum  detection  system  by 
applying  the  optimum  decision  rule  design  criterion.  In  fact,  one  would 
probably  be  satisfied  with  the  best  workable  system  of  several  possible 
choices.  What  is  needed,  then,  is  some  measure  of  individual  system  per¬ 
formance  which  may  be  used  to  compare  different  systems  and  which  will 
permit,  in  practice,  the  selection  of  the  best  system.  Such  a  performance 
measure  should  be  derivable  from  the  Neyman-Pearson  criterion  or  the 
essentially  equivalent  single  signal  threshold  criterion.  One  usable 
performance  criterion  is  discussed  in  the  following  section. 

In  summary  of  the  previous  subsections,  a  general  performance 
criteria  has  been  derived  in  the  form  of  a  decision  rule  which  minimizes 
the  expected  cost.  This  criterion  is  shown  to  be  of  limited  usefulness 
since  it  is  required  to  know  the  a  priori  probability  of  a  target  occurrence 
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and  cost  functions  assigned  to  the  various  outcomes.  The  Neyman -Pears on 
criterion  is  shown  to  be  free  of  these  limitations  and  yet  to  retain 
much  of  the  rigor  of  the  more  general  criterion.  Finally,  it  is  argued 
that  in  most  cases  the  Neyman-Pearson  criterion  corresponds  to  the  simple 
decision  criterion  of  setting  a  single  signal  threshold. 

e . .  A  Usuable  Comparison  Criterion 

It  is  desired  to  develop  an  analytical  measure  of  detection 
system  performance  which  may  be  used  to  compare  different  systems.  The 
various  systems  will  differ  in  their  decision  rules  since  the  decision 
rule  used  reflects  the  type  of  discrimination  scheme  of  a  particular  system. 

Methods  for  achieving  background  discrimination  are  briefly: 

(1)  Spatial  discrimination  by  using  small  instantaneous  fields 
of  view.  Thus,  if  the  resolution  element  (an  individual 
detector  of  an  array)  is  on  the  order  of  the  target  size, 
the  field  of  extensive  backgrounds  is  reduced.  This 
method  has  been  considered  in  the  single  scanning  detector. 

(2)  Spatial  discrimination  by  using  inherent  background  patterns 
and  correlations.  The  fact  that  the  spatial  frequencies 

of  cloud  radiance  differ  from  that  of  the  target  (essentially 
a  point  source)  might  be  utilized  for  discrimination.  Also, 
correlations  due  to  regularities  in  cloud  size  and  distri¬ 
bution  may  permit  discrimination.  An  example  has  been 
considered  for  clouds  several  resolution  elements  wide  in 
which  signals  from  adjacent  elements  are  compared  and  no 
alarm  is  sounded  if  both  signals  exceed  the  threshold. 

(3)  Spectral  discrimination  by  choosing  the  most  appropriate 
spectral  region.  It  appears  that  the  4.3  micron  band 
offers  no  great  advantage  over  the  2,7  micron  region,  how¬ 
ever. 
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(4)  Spectral  discrimination  by  examining  spectral  line  shapes 
or  by  comparing  fluctuations  in  two  or  more  different 
spectral  regions. 

(5)  Temporal  discrimination  by  determining  the  time  history 
of  a  potential  alarm  radiance.  Target  motion  or  self 
modulation  may  reveal  its  true  identity. 

(6)  Discrimination  accomplished  by  limiting  the  chances  of 
observing  a  severe  background.  Thus,  a  reduction  in 
coverage  say  by  limiting  the  scan  to  large  angles 
relative  to  the  sun  position  to  reduce  seeing  bad  for¬ 
ward  scattering  or  limiting  the  scan  to  regions  near 

and  above  the  horizon  to  take  advantage  of  transmissivity 
losses  might  be  useful. 

(7)  Polarization  discrimination  by  recognizing  that  radiation 
reflected  from  clouds  may  be  polarized  differently  than 
target  emission. 

A  particular  system  may  perform  measurements  upon  both  target 
and  background  signals.  In  general,  a  system  may  be  required  to: 

(1)  Detect  the  presence  of  a  completely  known  signal  in  the 
presence  of  background  noise. 

(2)  Detect  the  presence  of  an  incompletely  known  signal. 

(3)  Detect  an  incompletely  known  signal  and  measure  its 
unknown  characteristics. 

As  a  result,  the  decision  rule  may  become  an  involved  function 
of  various  system  operations. 

The  gain  functions  associated  with  the  correct  decisions  (C^  and  C^) 
will  be  neglected  and  the  system  which  produces  fewest  errors  of  target  miss 
variety  for  a  given  rate  of  false  alarm  error  will  be  considered  the  best 
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system , which  is  in  accordance  with  the  Neyman-Pearson  criterion.  It  is 
assumed  that  a  tolerable  false  alarm  rate,  represented  by  may  be 
assigned.  For  comparison  purposes,  systems  will  be  examined  when  their 
false  alarm  rates  are  identical.  Then  the  system  with  the  smallest  condi¬ 
tional  probability  of  miss  is  the  correct  choice.  Since 

^ra  “  / |t  0  (^1  -  5(S,n  )}  dS  -  1  -  (4.26) 

this  is  the  same  decision  rule  for  the  best  system  as  selecting  the 
system  with  the  largest  P^  from  all  systems  having  equal  P^. 

The  particular  system  which  is  best  for  one  value  of  P^  may 
not  be  the  best  at  all  other  values  of  P^.  Thus,  curves  of  P^  versus 
P^  for  various  systems  may  cross  at  some  point, and  the  order  of  best  per¬ 
formance  may  become  inverted. 

It  is  necessary  at  this  point  to  consider  how  the  s^  signals 
from  the  n  sensors  are  combined  to  yield  the  total  signal  S.  The  decision 
operator  will  be  determined  by  this  operation.  For  example,  all  s^ 
signals  might  be  summed  and  then  passed  through  a  signal  threshold  device, 
or  each  s  could  be  individually  passed  through  a  threshold  device  before 
combination.  Various  relationships  between  the  individual  signals 
might  be  required  as  a  condition  for  decision  making.  One  such  relation¬ 
ship  could  be  a  statistical  correlation  and  certain  restrictions  upon 
the  degree  of  correlation  could  make  up  a  decision  rule.  The  procedure 
whereby  the  s^  signals  are  processed  to  give  the  resultant  "S'  which  is 
involved  in  the  decision  operator  will  be  critical  and  will  be  the  feature 
of  the  various  systems  which  may  be  examined. 

The  decision  operator  must  be  modified  to  permit  the  inclusion 
of  the  procedure  for  combining  the  signals  which,  in  turn,  depends  upon 
the  discrimination  scheme.  Call  this  decision  operator  5(^(®£)>  • 
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Certain  simplifications  may  be  possible  in  the  expansions  of 
the  conditional  probabilities  pC^S  |  T  »  0)  and  p(S  |t  k  0)  appearing  in 
the  general  equations  for  and  (Eqs.  (4.19)  and  (4.20)).  For  ex¬ 
ample,  false  alarms  result  only  when  an  alarm  is  given  but  no  target  is 
present.  However,  a  target  may  be  present  in  the  i^^  signal  and  be 
missed,  yet  an  alarm  signal  sounded  due  to  some  other  s^  satisfying  the 
alarm  condition  even  though  s^  contains  no  target  signal.  This  fortuitous 
alarm  will  be  an  extremely  rare  event  since  it  requires  both  that  a  target 
be  missed,  which  for  most  systems  is  very  unlikely,  together  with  a 
false  alarm  event,  which  is  also  rare.  Thus  terms  in  p('s  |  T  “  0)  of  the 
sort  p(s.  i  W  .  >  s  .c  w  .  t.  ^0,  t .  =0)  may  be  neglected.  Here 
^  is  the  subspace  of  all  s^  for  which  a  decision  for  target  presence 
is  made. 

An  alarm  should  be  given  when  one  or  more  targets  appear. 

Because  of  the  spacing  of  missile  launches  in  both  time  and  space,  it  is 
unlikely  that  more  than  one  target  will  be  detected  at  any  one  instant 
although  the  probability  does  increase  as  the  spatial  separation  of  the 
i^^  and  sensors  becomes  larger.  It  is  probably  quite  unlikely  that 
two  missiles  will  be  launched  simultaneously  side  by  side  so  that  target 
signals  appear  in  adjacent  sensor  signals.  Thus,  terms  in  the  expansion 
of  p(S  I  T  ^  0)  of  the  sort  p(s^cQ  ^  i+1  ^  ^  ‘^i+l  ^ 

be  neglected. 

The  decision  rule  will  determine  the  detailed  form  of  Eqs.  (4.19) 
and  (4.20).  For  example,  suppose  a  system  of  two  sensors,  s^^  and  S2  and 
a  decision  rule 


6(s,  Q  ) 


^  1  when  >  Sg 

/  1  when  >  Sg 

0  when  ^  s^;  4  Sg 


(4.27) 
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Then 


(4.28) 


and  T  »  0  means  t|^  »  0  and  t2  *=  0.  No  assumption  about  the  correlation 
between  s^  and  S2  need  be  made  since  the  result  may  be  written  for  the 
general  case.  Then  is 


Pf  *  P(Sj^>S2  I  t^  *  0,  t2  =  0)  ^(s^,S2;  Sq)  ds^  ds^ 


(4.29) 


p(Si  |tj^  =0)  p(s2  I  t2  =  0,  s^)  dSj^  ds2 


'0 


Similarly  we  can  find  the  detailed  expression  for  P^.  T  ^  0  means  that 
(tj^  ¥  0,  t2  =  0);  (tj^  =0,  t2  ¥  0),  or  (tj^  ¥  0,  t2  ¥  0). 

Thus,  P^  becomes 

Pj  =  JJl^ p(Sj^,  S2)  I  tj^  V  0,  t2  =  0)  + 

P(Si,  S2  1  =  0,  t2  ¥  0)  + 


(4.30) 


where 


and 


P(Sl»  ^2  I  ¥  0,  ¥  oj  S(Sj^,S2;Sq)  dSj^ 

P(s. ,  Sj  I  t.  ¥  0,  tj  «=  0)- 

p(s.  I  t.  ¥  0)  p(Sj  |t.  ¥  0,  tj  =  0,  s.) 

p(Si,Sj  I  ti  ¥  0,  tj  ¥  0)  = 

P(s^  I  ¥  0)  p(Sj  I  ¥  0,  tj  ¥  0,  s^) 


ds. 


(4.31) 


(4.32) 
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Eqs .  (4.29)  and  (4.30)  above  take  into  account  correlations  between  back¬ 
ground  signals,  b^  and  b^,  and  correlations  such  as  the  probability  of  a 
double  target  (t^  ^0,  t^  ^0). 

Finally,  it  is  necessary  to  specify  how  b^  and  t^  components 
add  to  yield  the  resultant  s. .  If 

s.  -  b.  +  t.  (4.33) 

111 


and 

p(s. /t.  =  0)  «  w  (s.) 
^'ii  '  o'l' 


(4.34) 


then 


P(Si/tj^  k  0)  = 


(bp 


(s, 


b.)  db. 

X  1 


(4.35) 


/ 


w.  (t . )  w  (s .  -  t . )  dt , 
110  1  1  1 


Under  certain  conditions,  the  convolution  integrals  above  may  be  simplified 
by  approximation  if  some  of  the  properties  of  ^^i^  known. 

For  example,  if  the  probability  density  function  of  background  radiances 
is  less  than  some  small  fraction  of  the  comparable  missile  radiance  densities, 
then  a  good  approximation  is 


(sp 


(s  .  -  b. )  db.  CS  w,  (s  .  )  . 

'  1  i'  1  1  1^ 


(4.36) 


It  is  permissible  to  use  this  approximation  for  most  of  the 
background  distributions  considered  here.  Note  finally  that  the  background 
must  be  represented  as  an  n-dimensional  joint  probability  density  function 
given  by  Eq.  (4.37) 


(B)  =  p  (bp  p  (b^  I  bp  p  (b3  I  b^,  bp  P  (b^  I  b^.  b^ . ,b^_p 

It  is  seen,  then,  that  the  approach  used  here  determines  how  the  background 
is  to  be  described.  Unfortunately,  for  correlated  backgrounds  with  n  any 
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larger  than  2  or  perhaps  3  or  4,  the  practical  measurement  and  determination 
of  such  joint  probability  density  functions  becomes  well  nigh  impossible. 

The  following  section  considers  this  problem  in  more  detail  and  typical 
probability  density  functions  for  non-Gaussian  distributions  are  presented. 

4.3  PROBABILITY  DENSITY  FUNCTIONS 

a.  Properties  of  Correlated  Joint  Density  Functions 

In  Section  4.2  a  the  nomenclature  for  n~dimensional  probability 
density  functions  was  introduced  for  background  signals,  target  signals, 
and  resultant  signals.  It  is  probably  not  possible  to  obtain  these  complete 
joint  density  functions  for  n  any  larger  than  3  or  4.  If  the  n  random 
variables  were  statistically  independent,  the  joint  density  function  could 
then  be  written  down  as 


p(S)  =  p(sp  pCsp  ....  p(s^^) 


(4.38) 


However,  these  individual  random  variables  are  not  statistically  independent 
and  Eq.  (4.38)  cannot  be  used  to  represent  the  joint  probability  density 
with  good  accuracy.  Eq.  (4.37)  might  be  used  if  corrections  for  random 
variable  correlation  could  be  incorporated. 


Suppose  a  joint  probability  density  function  exists  for  two 

correlated  random  variables  s.  and  s.  designated  by  p(s. ,  s  .) •  Let  s. 

1  j  ®  ^  ^  1  y  1 

and  s,  have  individual  probability  density  functions  p(s,)  and  p(s.), 
respectively.  Also,  let  p(s^  I  j^  ^^^j  I ^i^  corresponding  condi¬ 
tional  probability  density  functions  for  s^  given  Sj  and  s^  given  s., 
respectively.  Then 


p(s..  Sj) 


p(Sj  I  s^) 


(4.39) 
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Suppose  next  that  the  correlation  between  and  is  perfect 
so  that  the  correlation  coefficient  P  is  unity.  Then  p(s^,  s^) 
becomes 


P(Si,  Sj)  =  p(s.)  5(s,  -  S^)  =  p(Sj)  5(Sj-S.) 


(4.40) 


where 


(5  (s.  -  s.)  and  u(s.  -  s.)  are  delta  functions 
1  J  ^  j 


5^(3.  -  s.)  «  ^(s .  -  s.)  = 

1  j  J  ^ 


1  s.  =  s  ^ 

^  j 


0  s.  V  s  . 
2.  ^  J 


(4.41) 


Finally  suppose  the  two  random  variables  are  statistically 
independent  so  ^  =  0,  Then  p(s^j  Sj)  is  given  by  Eq.  (4,42)  as 


P(s^>  Sj)  =  p(s.)  p(Sj) 


(4.42) 


By  comparison  of  Eqs.  (4.41)  and  (4.42)  with  the  general  form 

of  Eq.  (4.39),  it  will  be  seen  that  the  conditional  probabilities, 

p(s^  1  Sj)  and  p(Sj  j >  range  in  functional  structure  from  delta  functions 

for  perfect  correlation  to  the  individual  density  functions  for  no  correlation. 

Thus,  for  some  intermediate  value  of  P  ,  the  delta  function  in  Eq.  (4.40) 

must  be  replaced  by  a  function  broader  in  (s.  -  s.)  whose  dependence  upon 

1  J  p 

the  conditional  variable,  say  in  this  case  s^,  becomes  less  as  /  approaches 
zero.  The  actual  function  will,  in  most  cases,  be  difficult  to  estimate, 
but  the  above  considerations  indicate  something  of  its  general  form  as  a 
function  of  the  degree  of  correlation. 

The  n-dimensional  joint  probability  density  function,  p(S), 
must  satisfy  certain  properties  such  as 


p(S)  dS  =  1 


(4.43) 
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where 


p(S)  ds  .  .  .  .  ds  =  p(s  ) 
i  n-i  n 


p(S)  dS  =fjL 


p(S)  dS  =  +  ij^ 


s ,  s .  p (S)  dS  =  ^  ,  +  / i 

i  J  ij  M. 


(4.44) 


(4.45) 


(4.46) 


(4.47) 


(4.48) 


^  ■  M  )  -/W  )  > 


(4.49) 


are  expectation  values.  Note  that  p(S)  defined  by  Eq.  (4.38)  satisfies 
all  these  properties  except  Eq.  (4.49)  unless  ^  that  is,  unless 
there  is  no  correlation. 

It  may  in  some  instances  be  possible  to  approximate  p(S)  when 
correlation  is  involved  by 


p(S)  =  p(Sj^)  P(S2)  .  .  .  p(s^)  +  f(S) 


(4.50) 


where  f(S)  must  satisfy  the  properties 


f  (S)  dS  =  0 


f(S)  ds.  .  .  .  ds  1=0 
^  ^  1  n-1 


s^  f(S)  ds  *  0 


s^^  f(S)  dS  =  0 


(4.51) 


(4.52) 


(4.53) 


(4.54) 
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s,  s ,  f(S)  ds  *  OC  ,  , 

1  J  ij 


(4.55) 


The  problem  reduces  to  finding  the  appropriate  form  for  f(S).  In  general, 
the  approximation  for  p(S)  of  Eq.  (4.50)  will  not  satisfy  higher  order 
properties  of  the  true  p(S). 

In  summary,  two  approaches  to  approximating  joint  probability 
density  functions  in  the  case  of  correlation  have  been  suggested,  namely, 
coLj  -ct  successive  P(Sj^)  in  Eq.  (4.38)  to  account  for  the  correlation 
or  assume  independence  and  then  make  an  additive  correction  (Eq.  4.50). 

b.  Deduction  of  Probability  Density  Function  from  Physical  Models 

If  physical  models  for  the  background  radiance  and  target  radiance 
can  be  derived,  say  from  theoretical  and  empirical  studies,  how  may  the 
required  probability  density  functions  be  obtained?  In  general,  a  set  of 
particular  conditions  must  hold  for  the  production  of  a  given  radiance 
value.  If  p^  is  the  probability  oi  the  occurrence  of  the  i^^  condition, 
and  if  the  conditions  are  statistically  independent,  the  probability  of 
obtaining  the  new  radiance  is 

p(K)-^  TTiPik 

k 

where  the  sum  of  k  is  over  all  sets  of  i  which  separately  produce  H. 
Although  H  is  a  function  of  all  i,  each  p^  will  be  determined  by  the 
specification  of  H  and  the  particular  set  k  giving  rise  to  H .  Finally, 
if  the  p^  are  not  statistically  independent,  a  convolution  of  probabilities 
will  be  required. 

In  actual  practice,  p(H)  will  probably  have  to  be  obtained  from 
the  physical  model  by  Monte  Carlo  computer  technique.  It  may  also  be 
easier  to  obtain  an  estimated  cumulative  distribution  function, P  (H)  and 
then  differentiate  with  respect  to  H  to  obtain  p(H). 


4-27 


At  this  point  another  difficulty  arises  in  the  practical  applica¬ 
tion  of  the  method  for  assessing  system  performance.  It  has  been  shown 
how  the  probability  density  functions  for  model  cloud  background  radiances 
may  be  obtained  in  principle.  However,  Eqs.  (4.29)  and  (4.30)  show  that 
the  probability  density  function  of  the  signal  as  it  appears  to  the 
decision  operation  is  what  is  actually  required.  The  cloud  radiance  is 
converted  by  the  sensor  to  a  voltage.  At  the  same  time  some  modification 
of  the  original  radiance  probability  density  function  may  occur  since  the 
sensor  may  not  be  truly  linear  for  all  radiance  values  or  previous  signals 
may  have  an  effect.  Also,  the  sensor  may  introduce  noise  of  its  own  which 
definitely  changes  the  density  function.  As  the  signal  processing  proceeds, 
the  linearity,  noise,  and  bandwidths  of  preamplifiers  will  further  change 
the  signal  distribution.  The  signal  presented  to  the  decision  operation 
device  will  represent  possibly  a  much  different  probability  density  function. 
Unless  these  modifications  of  the  original  signal  are  well  understood, 
it  will  not  be  possible  to  relate  the  background  densities  to  system  per¬ 
formance.  We  will  show  immediately  below  what  simplifying  assumptions  may 
be  made  about  the  system  to  remove  this  difficulty. 

Let  the  input  irradiance  to  the  i^^  sensor  be  H^(0,  0).  Because 
of  the  scanning  law  and  sensor  motion 


0  -  0  (t) 
0  “  0  (t) 


(4.57) 


so  that 


H^(e,  0)  =  H.(t) 


(4.58) 


It  is  assumed  that  the  probability  density  function  P(H^)  known  and 
that  it  is  required  to  obtain  p(s^). 

The  output  signal  at  time  t  corresponding  to  H^(t)  is  s^(t). 
However,  contributions  to  s^(t)  result  not  only  from  the  instantaneous 
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value  of  H^(t)  but  also  from  the  past  values  of  H^(t).  This  is  due  to 
the  system  response  function  and  its  inherent  rise  time  and  causes  the 
difficulty  in  converting  p(s^)*  Another  source  of  difficulty 

is  the  noise  introduced  by  the  system.  However,  we  assume  that  the 
internal  system  noise  is  sufficiently  low  so  as  to  be  neglected  without 
appreciable  error. 

The  way  in  which  the  history  of  is  taken  into  account  is  to 
de.i.ie  a  system  impulse  response  function,  k(t) ,  which  when  convoluted 
with  H.  gives  s  .  ,  i . e . 

‘  r 

s  (t)  =  i  H  (t  -  A  )  k  (y()  d  X  •  (4.59) 

However,  this  simple  integral  gives  the  correct  signal  only  for  linear 
systems.  Higher  order  integral  equations  are  required  for  other  types 
of  signal  processing.  In  our  case,  we  are  interested  only  in  the  signal 
processing  which  occurs  up  to  the  input  to  the  decision  operator,  and 
this  processing  is,  in  many  cases,  closely  linear. 

Eq.  (4.59)  relates  individual  values  of  H.  occurring  over  a 

t  h  ^ 

sampling  time  interval  of  the  i  sensor  to  the  output  signal,  s^. 

However,  what  we  need  is  some  relationship  between  pCs^)  and  p(H^),  not 
between  individual  values  of  the  random  variables.  In  general,  such 
a  relationship  cannot  be  found  analytically  except  for  certain  special 
or  simple  cases  such  as  Gaussian  background  distributions  or  for  sensors 
with  sufficiently  rapid  rise  times. 

Suppose  the  impulse  response  function  delays  the  signal  by  an 
amount  t^  but  otherwise  doesn’t  change  the  signal.  Such  a  function  is 
represented  by  the  delta  function 
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For  this  case 


s.(t)  -  H.  (t  ) 
1  L  o' 


and  the  probability  density  functions  are  identical,  i.e.. 


(4.60) 


P(H^)  =  P(s^) 


(4.61) 


Finally,  if  is  a  slowly  varying  function  of  time  relative 

to  ^he  time  constant  of  the  detection  system,  i.e.,  ^oes  not  change 

appreciably  over  a  time  interval  from  t  -  7^  to  t  +  7^  where  'f  is  the 
time  constant  of  the  system,  then 


s.  (t) 


H.  (t) 


fOO 


k  ( A  )  dA 


(4.62) 


The  problem  of  obtaining  P(s^)  vanishes  for  this  special  case.  Actually, 
all  that  is  required  is  that  the  relationship  between  s^(t)  and  H^(t) 
be  raonotonic  so  that  a  single  value  of  s^(t)  exists  for  a  given  H^(t). 
Then,  again,  Eq.  (4.61)  is  valid. 

For  a  high  frequency  cutoff  filter,  the  transfer  function  is 


-rruJf 


(4.63) 


where  'T  is  the  characteristic  response  time  of  the  system.  The  Fourier 
transform 


k(t)  I  K  (") 


TT 

of  this  function  generates  the  impulse  response  function 

-t/r 


(4.64) 


Ke 


k(t)  = 


T 

0 


t  >  0 
t  <  0 


(4.65) 
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Let  us  Fourier  analyze  what  happens  to  the  signal  as  it  is 
processed.  The  Fourier  transform  of  the  input  signal  is 


.00 


H(t)  =  IT  I 

(°o 

%  (a/)  = 


-  \  ijJ  t 

H(t)e  J 


/ 


Eq.  (4.59)  becomes  upon  substitution 

foo 

Hit-  /V ) 

to 


S  .  (t)  =  ^ 

1  TV 


H(t-  diud^ 


(4.66) 


(4.67) 


/OO 


1 

r 


K.  i^)dij^  / 

/o  /  o 


H(t-  A)ej‘^^dA 


/OO 


(4.68) 


=  -^  I  /((a;)ej‘^‘'  dw  j  H(t-  d  A 


/OO 


1 

TT 


/<  (^)  M'(a))  duJ 


For  K.  (6U)  given  by  Eq.  (4.63)  the  product  {00)  l^(cu)  which  is 
effectively  the  Fourier  transform  of  s^(t)  will  be  small  for  high  fre¬ 
quencies  (large  values  of  cU  )  .  To  see  this,  expand  (4^  {Uj) 

2 


(tJ)  =  M  (0)  +  a;h  (0)  +  ^  H.  (0) 


(4.69) 


a) _ n 


(0) 
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For  large  CO  ,  multiplication  by  (U>)  effectively  reduces  the  power 
of  caJ  in  each  term  and  the  new  series  converges  more  rapidly  as  ^ 
increases  in  size  than  the  original. 

If  the  loss  of  high  frequency  components  of  H^(t)  is  small 
compared  to  the  mean  value  of  H^(t),  then  the  radiance  probability  density 
function  will  closely  approximate  the  signal  density  function.  If,  how¬ 
ever,  the  high  frequency  losses  are  high,  the  probability  density  function 
will,  be  appreciably  altered. 

If  k(t)  is  not  an  extremely  narrow  function  of  t,  but,  on 
the  other  hand,  is  not  excessively  broad  either,  then  an  approximation 
may  be  used  with  good  results.  Thus  for  each  s^(t)  given  by  Eq.  (4.59) 
there  will  correspond  an  interval  of  input  signals  (H^, 

The  width  of  this  interval,  depends  upon  the  width  of  the  weighting 

function,  k(t)  .  An  average  probability  density,  i-S  assigned  over 

this  interval.  Then,  approximately 

p(s.)  ds.  ^  p  S  H. 

1  ^ave  ^  L  (4.70) 

The  wider  k(t)  the  poorer  will  be  the  approximation  of  Eq.  (4.70). 

Another  point  of  view  with  regard  to  the  problem  of  relating 
input  and  output  probability  density  functions  is  simply  that,  when  the 
input  is  a  random  variable  rather  than  a  time  variable  of  known  analytic 
form,  no  Fourier  transform  exists  for  such  a  variable.  As  a  consequence, 
the  modification  of  the  input  density  by  the  system  cannot  be  evaluated 
unless  the  modification  is  negligible  or  slight  so  that  the  transformation 
of  the  random  variable  reduces  to  a  simple  algebraic  operation.  In  order 
to  handle  this  problem  when  the  optimum  linear  filter  is  considered  in 
Section  4.6,  the  random  functions  are  converted  to  autocorrelation  functions 
(with  consequent  loss  of  phase  information)  for  which  transformations  are 
possible. 
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c .  Typical  Background  Distributions 

Next,  consider  some  typical  examples  of  non-Gaussian  background 
distributions.  Such  distributions  are  required  in  order  to  account  for 
radiance  from  several  different  sources.  For  example,  partial  cloud 
cover  gives  rise  to  a  bimodal  distribution  of  radiance;  one  mode  is  cen¬ 
tered  around  the  clear  sky  mean  and  the  other  around  the  cloud  mean. 

Clouds  at  two  different  levels  may  be  responsible  for  a  bimodal  distri- 
bull v)n  because  of  the  differing  transmissivity  of  the  paths  to  and  from 
the  clouds. 

Some  properties  of  bimodal  noise,  each  mode  being  represented 
by  a  Gaussian  function,  are  examined  below.  The  underlying  assumptions  are 

(1)  Noise  and  signal  are  uncorrelated  from  one  resolution 
element  to  the  next. 

(2)  Noise  and  signal  are  additive. 

(3)  A  decision  is  made  that  signal  is  present  if  t  +  b  ^ 
or  absent  if  t  +  b  ^  s  where  t  and  b  are  signal 
and  noise  voltages  at  the  output  of  the  radiation 
detector  and  s^  is  a  threshold  voltage. 

(4)  The  noise  probability  density  function  is 

p(b)db  =1  yppf^  ~  cr^ 

^  (4.71) 

where  the  a^  and  2^  are  called  weights  of  the  distribution.  They  might 
physically  represent  the  probabilities  of  the  presence  of  high  or  low 
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clouds.  and  jji  ^  called  modes  of  the  distribution  and  ^2-^1 
is  called  the  spacing  of  the  distribution.  Physically  Ml  and  u  2  might 
represent  the  mean  radiance  of  high  and  low  clouds  multiplied  by  the 
responsivity  of  the  sensor,  the  size  of  the  aperture  stop  and  the  subtend 
of  the  field  stop,  and  called  the  variances  of  the  distri¬ 

bution.  Physically  they  might  represent  the  variance  in  radiance  of  high 
and  low  clouds  multiplied  by  the  same  factors  as  jji  and  IX2' 

The  false  alarm  probability  and  signal  detection  probability 
have  been  evaluated  under  the  five  assumptions  by  integrating  numerically  the 
noise  density  function  and  the  signal  plus  noise  density  function.  Some 
results  for  the  false  alarm  probability  are  shown  in  Figures  4.1  through  4.5. 
These  figures  show  false  alarm  probability  plotted  horizontally  from  right 
to  left  against  threshold  voltage  plotted  vertically.  The  graphs  are  on 
"probability  paper"  on  which  a  Gaussian  distribution  plots  as  a  straight  line. 

Figure  4.1  shows  the  effect  of  varying  the  weights  of  the  distri¬ 
bution  from  =0.5  to  a2  =0.01.  Figure  4.6  shows  a  graph  of  the  proba¬ 
bility  density  function  for  the  extreme  case  a^^  =  0.99,  =  0.01. 

Figure  4.1  also  shows  the  false  alarm  probability  characteristics  for  equiva¬ 
lent  Gaussian  distributions  --  that  is,  the  false  alarm  characteristics 
that  would  be  anticipated  by  a  system  designer  who  had  measured  the  mean 
and  variance  of  the  background  noise  and  had  fitted  it  with  a  Gaussian 
distribution.  These  characteristics  are,  of  course,  straight  lines. 

The  signal  detection  probability  characteristics  under  these 
simple  assumptions  are  obtained  by  displacing  the  false  alarm  curves  of 
Figures  4.1  through  4.5  vertically  by  the  signal  amplitude  in  the  same 
units  as  the  threshold  scale.  In  Figures  4.2  through  4.4,  the  weights  are 
a^^  =0.9  and  a^  =0.1. 

Figure  4.2  shows  the  effect  of  varying  the  spacing  of  the  two 
modes  by  holding  fJi  ^  fixed  and  progressively  decreasing  jJi2*  the 
spacing  decreases  the  characteristic  approaches  Gaussian  behavior. 
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FIGURE  4-1.  EFFECT  OF  VARYING  WEIGHTS 
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FIGURE  4-2.  EFFECT  OF  VARYING  SPACING 
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FIGURE  4-3.  EFFECT  OF  VARYING  WIDTH  -  a  VARIABLE 


iiiiiiiiiiiiiniiiBkiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!i!!! 

«SKrJKSSSSBi£»fSr.SSSSS!SS!!!!SSiS^^^^^^^ 

■SbSSSSSSS  SSSSSSS^SS  SSa S5SSSSSSS  8ggs:as:sg8;g:a8;as:8aa::g:85:g8;l 

5S5S555&555E585558B^V8HsviiBBMiSMS  viMM&MivSBtMB  mm  **22288228 


~^ss»s8Si:KS8ss:!::ss::BS8Ss»8»:s:f - * 


■■■■S55a««  ■■■■■■aaaaaaaaBaaaaa  maaaaHaaa  ■■■■■■■■■■  ■■■■■■■aaiMaaaaiiaaa 

■■■■■■■■MM  ■■■aaaaBmiBaaaaaHaaMMaaaMa'iaa 

■■■■■■■■■a  BaaBa ■■■■■ ■■■«■ aaBaa  BBaaaBBaaa aBaflBaaaaa  aaaaaBBBaa  aaaaa aar aa 

iiiiiiiiiiiiiiiiiiiiimiiiiiiiiniiiiiir - 

lllllllllllllllllllllillllllllllll 


I  SaaSSaSSaaSSaHiiBSaBSSSi 

I BBBBBBaBBBaaBBi - 

I BBBBiBBBaaBBaBi 


KBBBiBtfSSSSSa 

aaBBaaliBaSa 


I  ■HHgaH  WHiHH  HHHpgBBBBaaBHBBBBBtgim 

I  '■■■■■■■■■■  BMBBaBHBI  ■■■■■■■■■■■■■■ 


S3a 


iiUHiiiiu  iiiiiimiina 


ssKsSssSsIsSaiiSSai 

I  anHaKHMiaipB  ■■■«■  ■■■«■ 

I aBBBBBBBBB  BBBBBBMBBB _ 

laBBBBaaBBBaBBBaBBaBiBBMiaBM 

iiiiiiuiiiiiiiiiiiiiiiiiir 

iiiiiiiiiiiiiiiiHiiiiiiiii 


null 

lino 

ggsgSv 


HI! 

101 


HI 


liinniilimi 

10000000 


liHBBiBBBiiHBHBBBBBHBBBBBHiHli 

iiiiiiiiiiiiiiiiiiiimiiiiiii 

lllllllllllllllllilllllllllil 


'DNIIiaS  QlOHSaHHI 


4-38 


FIGURE  4-4.  EFFECT  OF  VARYING  WIDTH  -  VARIABLE 
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FIGURE  4-5.  TRIMODAL  GAUSSIAN 
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FIGURE  4-6.  BI-GAUSSIAN  PROBABILITY  -  DENSITY  FUNCTION 


Figure  4.3  shows  the  effect  of  increasing  the  width  of  variances 
of  the  second  term.  Figure  4.4  shows  the  effect  of  increasing  the  variance 
of  the  first  term.  Figure  4.5  shows  the  false  alarm  characteristics  of 
five  cases  of  rather  closely  spaced  trimodal  noise. 

In  general,  it  is  seen  that  the  false  alarm  characteristic  of 
a  simple  ''independent  look'*  system  operating  against  bimodal  noise  can 
be  substantially  different  from  the  characteristic  expected  from  a  similar 
system  operating  against  noise  which  has  been  assumed  to  be  Gaussian  even 
in  the  case  that  the  second  peak  is  small  enough  that  it  might  well  be 
missed  in  a  background  measurement.  The  bimodal  characteristics  approach 
the  characteristic  of  the  high  mode  at  very  high  thresholds  and  the  low 
mode  at  very  low  thresholds.  Between  these  two  extremes  is  a  transition 
hump.  The  magnitude  of  the  hump  decreases  as  the  weight  of  either  mode 
is  decreased  or  as  the  spacing  of  the  modes  is  decreased* 

Figure  4.7  shows  another  multimodal  background  radiance  proba¬ 
bility  density  function.  The  low  radiance,  due  to  atmospheric  scattering, 
blackbody  radiation,  and  reflection  from  very  low  clouds,  is  represented 
by  a  delta  function  centered  at  or  near  the  origin.  The  medium  radiance, 
due  to  possible  cloud  cover  at  the  tropopause,  and  the  high  radiance,  due 
to  occasional  very  high  altitude  clouds,  are  represented  by  two  Gaussian 
functions.  This  distribution  differs  from  the  above  bimodal  cases  in 
that  the  Gaussian  functions  have  much  larger  variances  relative  to  the 
mean  and  in  that  the  bulk  of  the  background  radiance  is  contained  at  low 
radiance  levels  (the  delta  function).  The  two  Gaussian  functions  then 
represent  the  structure  of  the  high  radiance  tail  of  the  total  background 
probability  density  function.  This  is  the  important  region  of  the  total 
density  function  since  it  is  the  high  background  radiances  which  must  be 
distinguished  from  targets.  The  equation  plotted  in  Figure  4.7  is 
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Wo(B)  X  10 
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where 


W^(b) 


-(b. 


i  —  I  ■> 

'yU  TTo 


+  32  ^  2>^^^<T-2 

VT^cr^ 


a 


o 


<^2 


0.94 

0.05 

0.01 


“6  2 

=  10  watt/cm  ster 

-5  2 

<r  =  10  watt/cm  ster 

2 


(4.72) 


Such  a  background  radiance  distribution  may  be  related  to  a 
distribution  of  cloud  heights  since 


P(h) 


,  .  db  .  .  ^>^0  dt 

w  (b)  -TT  =  w  (b)  — s  — rr 
o^'^dh  o^^  t  dh 


(4.73) 


where  p(h)  is  the  probability  of  a  cloud  at  an  altitude  between  h  and 
h  +  dh,  f^  ==  1/cos  i  -f*  1/cos  e,  t  is  the  transmissivity,  and  dt/dh  is 
obtained  from  Table  3.2. 


Under  the  assumption  of  a  coplanar  case  (f^  =  const.),  a  maxi¬ 
mum  overall  radiance  of  2.6  x  10"^  watt/cm^  ster  (see  Section  2,  Fig.  2.5), 
(i  =  75^,  and  9  =  37°,  Fig.  4.8)  is  obtained  for  the  cloud  height  distri¬ 
bution.  The  two  peaks  shown  correspond  to  altitudes  of  16  km  and  28  km. 

The  probability  density  function,  p(h),  has  been  normalized  to  unity. 

The  low  altitude  distribution  is  not  shown  in  detail,  but  there  is  a 
95%  probability  of  having  no  cloud  above  12  km  altitude  (26  kft) . 
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FIGURE  4.8*  CLOUD  DISTRIBUTION 
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4.4  SINGLE  SENSOR 

Consider,  now,  a  detection  system  consisting  of  a  single 
detector  element  and  having  no  memory  or  correlator.  The  response  function 
will  be  assumed  to  be  narrow  so  that  for  all  practical  purposes  the  proba¬ 
bility  density  functions  of  the  background  or  target  radiances  are  also 
good  approximations  of  and  Wj^(t),  respectively.  Let  the  decision 

for  alarm  be  the  condition  that  the  total  received  signal,  s,  exceed 
a  threshold  level,  s^.  In  this  simple  case,  is  given  by  Eq.  (4.19) 
in  which 

f 

S  =  (s)  , 

<  (4.74) 

T  =  (t) 

V 

T  =  0  means  t  »  0,  and 

>  Q  )  is  given  by  Eq.  (4.25)  so  that  Eq.  (4.19)  reduces  to 

p(s  I  t  =  0)  ds  (4.75) 

p(s  I  t  ^  0)  ds 
and  (4.36)  results  in 

w  (s)  ds 
o 

w^(s)  ds 


(4.76) 


(4.77) 

(4.78) 


Similarly,  is 


^d  = 


s 


e 


Using  Eqs.  (4.34),  (4.35) 

.00 


4 


"d  = 
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The  specification  of  P-  and  w  (s)  determines  s^  according 

to  y 

to  Eq.  (4.77).  This  value  together  with  w^(s)  then  determines  for 
the  case  of  the  single  sensor  (Eq.  (4.78)).  Figure  4.9  shows  P^  for 
w^(s)  given  by  Figure  4.7.  P^  may  be  obtained  for  any  desired  value 
of  if  the  target  probability  density  function  is  available. 

If  we  now  compare  the  result  of  our  performance  prediction 
(P^  for  a  given  P^)  with  the  prediction  which  would  have  resulted  from 
aii  ssumption  of  a  purely  Gaussian  radiance  distribution  matched  to 
the  actual  distribution  (or  simply  use  a  minimum  root-mean-square  noise 
criterion),  it  is  readily  apparent  that  the  structure  of  the  high 
radiance  tail  of  an  actual  distribution  has  a  pronounced  affect.  (See 
also  the  discussion  of  Section  4.3  c  and  Fig.  4.1).  The  structure  of 
the  high  radiance  tail  does  have  an  affect  upon  the  calculated  P^ 
which  is  essentially  ignored  under  the  r.m.s.  approach. 


If  w^(s)  and  Wj^(s)  are  known,  another  way  to  find  the  best 

Sq  for  the  single  sensor  is  to  choose  it  such  that  the  ratio  of  P./P^ 
y  Q  t 

obtained  from  Eqs.  (4.77)  and  (4.78)  is  a  maximum.  This  value  of 

s^,  ,  is  found  by  solving 

y  y  ,m 


w^(s)  ds 


w  (s)  ds 
o 


(4.79) 


or 


Wj^(s 

w  (s 
o' 


8,m 

0  ,m 


) 

) 


(4.80) 
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4.5  CORRELATED  SENSOR  OUTPUTS 

Substantial  gains  in  system  performance  are  achieved  when  the 
outputs  of  individual  sensors  are  correlated.  Suppose  a  detection  system 
consisting  of  two  adjacent  sensors.  The  decision  rule  is:  signals  from 
each  of  the  sensors  are  compared.  No  alarm  is  given  if  the  signals  from 
each  one  exceed  a  threshold.  An  alarm  is  sounded  if  the  signal  from 
either  one  exceeds  the  threshold  while  the  signal  from  the  adjacent  sensor 
do^  not.  Finally,  no  alarm  is  given  if  the  signals  from  each  sensor 
individually  are  below  the  threshold. 


A  false  alarm  may  well  be  the  result  of  two  erroneous  outcomes. 
The  signal  from  sensor  #1  may  exceed  the  threshold  while  #2  remains  below, 
or  vice  versa.  Since  the  cloud  image  will  almost  always  be  much  larger 
than  a  single  sensor  resolution  element,  it  is  unlikely  that  one  sensor 
will  view  an  infrared  reflecting  cloud  while  the  other  does  not.  Thus, 
the  probability  of  false  alarm  is  reduced.  Note  finally,  that  the 
occurrence  of  the  two  targets  in  adjacent  resolution  elements  is  assumed 
to  be  extremely  rare  and  so  is  neglected. 


The  decision  rule  is 

1  when  >  s^;  4  Sg 

1  when  >  Sg;  4  Sg 

0  when  >  Sg;  S2  >  Sg 

0  when  ^  Sg;  s^  <  Sg 

v 

The  assumptions  made  so  far  may  be  summarized  as  follows: 


5(s,  Q)  =  I 


(4.81) 
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r . 


S.(.,.  .b 


Sj  .  bj  + 


(Sj  I  t.  .  0)  -  hi 

<“l  I  ‘i  *  “>  ^ 


T  =  0  *  (t.  =  0,  =  0) 


T^O*  (t.  ^0,  t,  =0)  or  (t.  \  0,  t.  ^0) 
^  J  ^  J 


P  (t.  0,  t.  ^  0)  ^  0 


(4.82) 


Then 


p(S  I  T  =  0)  p  ^2  j  =  0,  =  0) 

“  P(Si|  =  0)  p(s2|  =0,  t2  =  0,  s^) 


(4.83) 


Also 


p(S  j T  ¥  0)  i  P(s^,  j  ^  0,  =  0)  (4.84) 

+  P(s^.  ^2  I  =  0,  ^  0)  +  p(s^,S2  j  %  0,  ^  0) 


But  it  is  assumed  tha^  the  condition  represented  by  the  last  term  on 
the  right  of  Eq.  (4.84)  is  impossible  so  that  this  term  must  be  excluded. 
Also,  because  of  symmetry 


p(S  |t  ^  0)  =2  p(s^,  I  k  0,  =  0) 

=  2  p(s^  I  ^0)  p(s2  I  0,  =  0,  s^) 


(4.85) 
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The  following  additional  assumptions  are  made: 

P(S2  j  ^  0,  =  0,  =  p(s2  j  =  0) 

i.e.,  P(s2  j  ^2  ~  independent  of  the  fact  that  t^  H  0, 

and 

P(S2  ^2  “ 

P(Sj^  I  ^  0)  =  Wj^  (Sj^) 

Using  the  above  information  and  Eqs.  (4,19)  and  (4,20)  gives,  finally 


/OO 


Pf  =  2 


"d  =  2 


/Sg  / 
/OO, 


WqCSi)  p(s2  jt^  -  0,  s^  =  b^)  ds^  ds^ 
82*0 


s  2  *“0 


/®0  / 

Eqs.  (4.88)  and  (4.89)  should  be  compared  with  Eqs.  (4.77)  and  (4.78) 
for  the  single  sensor.  The  modification  due  to  correlation  will  be 
apparent . 


The  ratio  of  P  ,/P^  is 
^  d  f 

too 


I  0 


Wq  (sp  ds^ 


OO  /S, 


0  '0 


Wq  (Sj^)  P  (s2  jt^  =  0,  =  bj^)  dSj^  ds. 


(4.86) 


(4.87) 


(4.88) 


(4.89) 


(4.90) 
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The  term  p(s2  j  t2  “0,  be  simplified  since  it  represents 

the  probability  of  finding  S2  knowing  that  the  adjacent  background  signal 
was  and,  therefore,  it  contains  the  background  correlation  information. 

A  comparison  of  Eq.  (4.90)  with  the  ratio  for  the  single 

sensor  discussed  in  the  previous  section  indicates  the  possibilities  for 
improvement  of  system  performance  using  correlation  techniques. 

a .  Improvement  due  to  Correlation 

It  should  be  readily  apparent  that  the  evaluation  of  Eq.  (4.90) 
for  any  but  exceptionally  simple  cases  would  be  extremely  difficult.  How¬ 
ever,  let  us  restate  Eq.  (4.90)  in  simplified  terms  and  introduce  some 
indicative  numbers  to  see  what  degree  of  improvement  is  possible.  For 
the  single  sensor 


and 


for 


single 


two  sensors  with 


p  (S  >  Sq  I  t  0) 
p  (s  >  ®e  I  t  =  0) 
correlat ion 


(4.91) 


P  (s^  ^  0)  P  1  ^  0,  =  0,  s^)  ^2) 

p  (Sj^  >  Sq  =  0)  P  (s^  4.  Sq  I  =  0,  =  0,  Sj^) 


The  improvement  is  represented  by  the  factor 


p  (^2  4  ^0  ^  0>  ^2  =  0,  Sj^) 

^  hi  "  *^2  “ 


(4.93) 


If  the  correlation  of  the  background  is  strong,  then  this  factor  will  be 
appreciably  greater  than  one.  Suppose,  for  example,  that  if  s^^  is  found 
to  exceed  s^  and  if  t^^  =0  and  t2  =  0  then  S2  has  a  probability  of  exceeding 
Sq  of  0.9.  Also,  suppose  that  if  s,  is  found  to  exceed  s^  but  t.  V  0  that 
the  probability  that  S2  will  exceed  given  t2  -  0  is  a  very  small  number, 

^  .  Then  the  improvement  is 
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I 


1  >  6 

1  -  0.9 


10 


An  example  of  the  relationship  between  a  simple  cloud  model  and  the 
improvement  will  be  examined.  Suppose  that  a  severe  background  is 
represented  by  a  cloud  10  detector  elements  wide.  We  then  ask,  what 
is  the  probability  that  if  the  i^^  detector  sees  the  cloud  the  i  +  1 
detector  will  not  see  the  cloud.  If  the  cloud  is  uniform  in  radiance, 
then  the  i^^  detector  must  be  viewing  near  an  edge  in  order  that  the 
i  4  1  element  not  view  the  cloud  also.  We  estimate  the  chances  of  this 
as  two  in  ten  or  a  probability  of  0.2.  In  this  case,  the  system  per¬ 
formance  would  be  increased  by  a  factor  of  5  because  of  correlation. 

A  cloud  20  elements  wide  would  give  the  improvement  factor  of  10 
obtained  above  for  the  conditions  assumed. 

b.  Correlation  Coefficient 

A  useful  measure  of  correlation  is  the  correlation  coefficient 
which  has  been  mentioned  before.  The  correlation  coefficient  for  the 
simple  cloud  model  introduced  in  the  preceding  section  is  derived  immedi¬ 
ately  below. 

The  correlation  coefficient  is  the  ratio  of  the  covariance  of 
two  random  variables  to  the  product  of  the  square  roots  of  their  variances. 


/ 


xy  = 


C  ov  (x.y) 


Vvar  (X)  Var  (y) 


(4.95) 


In  terms  of  the  joint  probability  density  function  of  x,  y, 
p  (x,y)  the  correlation  coefficient  is 


p  L 

/xy  =  — 


(x  -  M  )  (y-  U  '>  p  ("‘■y) 
X  y 


(4.96) 


(x-  p(x,y)  dxdy 
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What  is  the  correlation  coefficient  corresponding  to  the 
example  of  a  cloud  ten  elements  wide?  We  have 

E  (l>i  -  E  (b^  ^  j 

where  E  represents  the  expectation  operator.  We  take  for  our  model  a 
discrete  distribution  of  b. . 


P  (b.)  =  0.005 


E  (bp  -  E  ^  r.  P(bp 


(4.98) 

(4.99) 


and  assume  that  when  the  cloud  is  observed,  b.  =  10  which  exceeds  s 

1  0 

and  “  0  when  the  cloud  is  not  in  view.  Since  there  are  ten  positions 

where  this  value  of  b.will  be  received 

1 


now 


and 


E  (b.)  =  10  X  lO"^  X  .005  =  5  X  lO'^ 


(b.  -//)^  =  E  (b.^)  ^  b.^  ^  -  II 


/  9  -10  -12 

V  b,  P(b.)  =  10  X  10  X  .005  =  5  X  10 


1  '  1 


Also, 


(4.100) 


(4.102) 


E  (b.  -//  )  (b.  ^  ^-iU)  =  E  (b.  ^  ^)  -  /^  (4.103) 

and  since  there  are  nine  positions  for  both  b.  and  b.  ...  to  exceed  s^ 

^  L  1  +  1  0 


E  (b 


.  b.^p  °  ^  ^  X  .005  =  4.5  X  lO"^  (4.104) 
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Finally,  the  correlation  coefficient  is 


4.5  X  10~^^  -  2.5  X  10 

-12  -1 
5  X  10  -  2.5  X  10 


=  0.9 


(4.105) 


For  a  cloud  five  resolution  elements  wide,  a  similar  calculation  yields 
a  correlation  coefficient  of  0.8. 


4>  OPTIMUM  LINEAR  FILTER 

In  this  section  the  validity  of  using  optimum  linear  filters 
in  detection  systems  based  on  threshold  decision  criteria  is  examined. 
Some  simple  arguments  demonstrate  that  the  criteria  used  in  the  design 
of  optimum  linear  filters  do  not  usually  produce  optimum  system  per¬ 
formance  for  threshold  systems.  Also,  the  fact  that  the  exact  form  of 
the  target  signal  as  a  function  of  time  is  not  known  (the  target  signal 
iS  a  random  variable)  is  demonstrated  to  require  a  completely  statis¬ 
tical  approach  to  linear  filter  design. 

The  ideal  filter  would  be  one  with  no  attenuation  over  the 
frequency  band  occupied  by  the  target  signal  and  infinite  attenuation 
over  the  rest  of  the  frequency  spectra.  The  current  problem  is  to  decide 
on  the  best  signal  treatment  when  the  background  and  target  spectra 
overlap.  Internally  generated  noise  will  be  neglected  since  noise  due  to 
the  background  signal  levels  and  their  fluctuations  will  exceed  any 
internally  generated  noise. 

There  exist  several  possible  criteria  for  the  design  of  linear 
filters  such  as; 

(1)  Minimize  the  mean  square  error  at  the  output  of  the  filter 
(maximize  the  fidelity  of  target  signal  reproduction). 

(2)  Maximize  the  ratio  of  signal  power  to  average  background 
power  at  the  output. 
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For  the  purposes  of  target  detection,  the  latter  criterion  is 
used  since  this  would  be  expected  to  optimize  the  detection  sensitivity. 
When  it  is  required  to  retain  the  form  of  the  input  signal  with  the  best 
possible  fidelity  of  reproduction  in  the  output,  the  criterion  of  mini¬ 
mizing  the  mean  square  error  in  the  signal  at  the  output  is  applied.  The 
actual  output  signal  usually  consists  of  the  modified  input  signal 
together  with  some  passed  background  signal.  The  criterion  states  that 
the  mean  square  of  the  difference  between  the  actual  output  signal  and 
the  modified  ideal  signal  output  is  to  be  minimized.  The  ideal  output 
is  obtained  from  a  knowledge  of  the  pure  input  signal  or  its  spectral 
density  and  the  system  transfer  function. 


An  optimum  linear  filter  may  be  found  using  either  criterion 
regardless  of  the  background  probability  density  distribution  function. 
The  procedure  invariably  results  in  an  optimum  linear  filter  but  no 
assurance  is  given  that  a  more  optimum  non-linear  filter  does  not  exist 
Indeed,  for  non-Gaussian  distributions  of  background  with  weighted  high 
radiance  tails,  a  better  non-linear  filter  may  generally  be  found. 


Let  us  find  a  matched  linear  filter  for  a  non-Gaussian  back¬ 
ground  distribution  and  a  target  signal  of  known  time  dependence.  The 
matched  filter  is  one  where  the  spectral  density  (power  per  unit  band¬ 
width)  of  the  background  is  constant  for  all  wavelengths  so  that  the 
filter  is  "matched"  to  the  target  signal.  The  spectral  density  of  back¬ 
ground  is  not  generally  constant  so  the  matched  filter  must  be  preceded 
by  a  filter  which  converts  the  background  input  to  white  noise.  At 
the  same  time,  the  target  signal  input  will  become  modified.  Figure  4.10 
shows  the  hypothetical  system  being  considered.  The  input  background 
noise  is  b(t)  with  a  spectral  density  of  given  by 


(u>)  a 


)  dr 


(4.106) 


where  R,  (t)  is  the  autocorrelation  function  of  the  background. 
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FIGURE  4.10.  MATCHED  FILTER 


The  autocorrelation  function  of  the  background  input  signal  is 
required  in  order  to  obtain  the  spectral  density  of  the  background  (these 
two  functions  are  related  by  the  Wiener  theorem).  If  the  background 
irradiance  is  observed  for  a  period  of  tiine,T  the  autocorrelation  function 
is 

y  ^  HX+T)  d  (4.107) 

In  c^ctual  practice,  since  the  period  of  observation  does  not  approach 
infinity,  h(-T)  can  be  estimated  for  a  stationary  background  by 

/T  -  T 


Rb^  (T)  I  HX)HX+T)dX  (4.108) 

The  radiance  from  clouds  is  probably  best  described  as  a  random 
function  of  position  coordinations  with  some  provision  made  for  rough 
positional  correlations  which  may  exist.  For  such  a  case,  the  determi¬ 
nation  of  the  autocorrelation  function  becomes  a  statistical  problem.  The 
autocorrelation  function  is  defined  in  Eq.  (4.107)  as  a  time  average. 

There  exists  a  theorem  which  states  that  the  autocorrelation  function  may 
be  expressed  as  an  ensemble  average  (average  over  a  large  collection  of 
possible  background  patterns)  as  well. 

Because  of  the  practical  equivalence  of  the  time  average  and 
the  ensemble  average,  the  autocorrelation  function  can  also  be  expressed 


(4.109) 


Rb  (T)  - 


no®  fio,oo 

®1  ®2  ‘^^l  “  j  j  l®l)^ 


where  p(Sj^,  82  j"I“)  is  the  probability  that  Sj^  and  82  will  be  found  when 
the  time  eeparation  in  their  respective  measurements  is  T*  •  The  signals 
Sji  and  S2  may  be  considered  to  be  the  outputs  of  a  sensor  at  time  t^  and 
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respectively.  If  the  background  radiance  pattern  is  a  weak 
function  of  time,  and  may  also  be  considered  as  signals  from 
separate  identical  sensors  separated  in  the  object  space  by  the 
distance  ^  times  the  scanning  rate. 


The  average  background  power  out  of  the  prewhitening  filter 


will  be 


fOO 


h  (  TO 


J-  =  \  (tt))  W)  d 


cO 


(4.110) 


where  (<c?)  is  the  transfer  function  of  the  prewhitening  filter  and, 
in  general,  is  complex.  Kj^(cU)  is  its  complex  conjugate.  must  be 
chosen  so  that 


S,1  (O))  .  (4.1U) 

where  (a))  is  the  spectral  density  of  the  background  signal  at  the 
output  of  the  prewhitening  filter  and  is  a  constant  determined  by 
practical  limitations  of  real  sensors  and  preamplifiers. 

(a>)  =  i<0)  (CO)''  (cO)  =  (4.112) 

Thus,  with  a  knowledge  of  (a> )  and  the  choice  of  it  is  possible 

to  solve  for  the  (^ )  required  to  whiten  the  background. 

Since  the  autocorrelation  function  is  an  even  function,  i.e., 

h  (T)  =  Rb  (4.113) 

the  spectral  density  S^(^)  defined  by  Eq.  (4.106)  will  be  real,  and  it 
will  always  be  possible  to  find  an  appropriate  K^(<i>)  which  will  in 
general  be  complex. 
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The  input  target  signal,  t('X'),  will  also  have  been  modified 
by  the  prewhitening  filter  to  a  new  target  signal,  t^  {*C)  •  If  M  (W) 
is  defined  as 


r 

M  (UJ)  =  I 

/  -oa 


rjoJT  , 


(T)  dT* 


and 


t  (T)  = 


r 

TT  I 

/o 


then  it  may  be  shown  that 


ti  (r)  -:fjr 


Jb 


icur 


2  j  ^  M  (fO)  duJ 


M  (OJ)  (O))  d<o 


(4.114) 

(4.115) 

(4.116) 


The  criterion  of  maximum  signal  power  to  background  power 
ratio  requires  that  >  the  transfer  function  of  the  second  filter, 

be  chosen  to  make  a  maximum  of 


[  4  ] 

[  l>2  <T)  ] 


,oo 


(4.117) 


M  (u))  Kj^(W)  (O))  du) 


TT  S 


bl 


.1  -r. 

r 


(£l>)  K2  (fit))  doJ 

It  is  possible  to  show  that  the  appropriate  filter  transfer  function  is 


K2  (<0)  =  M*  {<*?).  (4.118) 

To  summarize,  what  has  been  done  above  is  to  apply  optimum 
linear  filter  theory  to  a  non-Gaussian  background  which  was  described 


4-59 


appropriately  for  this  approach  in  terms  of  the  autocorrelation  function. 
The  £ :tocorrelation  function,  spectral  density  (Wiener  spectrum),  and 
joint  probability  density  functions  are  related  background  descriptions 
and  their  use  in  deriving  the  optimum  linear  filter  was  illustrated. 

The  particular  approach  chosen  applied  the  criterion  of  maximum  signal 
to  background  signal  ratio. 

The  assumptions  which  were  made  were  that  the  target  signal 
W'  a  known  function  of  time  and  that  the  autocorrelation  function  of 
the  random  background  signal  was  given.  If  the  target  signal  is  also 
a  random  variable  then  a  matched  filter  cannot  be  obtained  by  appli¬ 
cation  of  Eq.  (4.118).  Missile  targets  will  generate  signals  with  a  random 
distribution  of  amplitudes  although  the  time  variation  of  the  signal 
may  actually  be  determined  by  the  scanning  mode  and  scan  rate  parameters. 

In  the  case  where  the  target  signal  is  a  random  variable  its 
autocorrelation  function  must  be  obtained  (Eq.  (4.109))  since  the 
actual  form  of  the  signal  is  not  known  with  certainty  but  only  pre- 
dictively.  If,  then,  the  target  Wiener  spectrum  (Fourier  transform  of 
the  target  autocorrelation  function)  overlaps  the  Wiener  spectrum  of 
the  background,  the  proper  choice  of  a  filter  is  not  immediately  obvious 
(Fig.  4.11).  However,  a  criterion  requiring  the  filter  to  generate  a 
maximum  root-mean-square  target  signal  to  root-mean-square  background 
signal  is  evidently  applicable.  Thus,  we  are  to  find  the  linear  filter 
(or  transfer  function,  K(Cc)))  which  maximizes  the  ratio  of  Eq.  (4.119). 


t 

out 

— Z 
b 

out 


s^(tJ)  K  (a>)  duJ 


(«J)  K  (cj)  dtJ 


where  (U))  is  the  target  spectral  density. 


(4.119) 
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FIGURE  4.11.  POWER  SPECTRA 
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Actually,  more  will  be  required  of  K(a^)  than  that  it  lead 
to  a  maximum  value  for  Eq.  (4.119).  For  example,  the  condition  of  a 
maximum  of  average  target  to  average  background  power  could  be  satis¬ 
fied  if  K(u))  is  infinite  over  that  range  of  where  is  zero 

but  S^(u))  remains  finite.  Such  a  filter  cannot  be  realized  in  practice 
since  an  infinite  filter  gain  is  implied.  It  should  be  possible  to 
find  some  compromise  K(U))  which  is  practically  realizable  and  yet 
permits  good  signal  processing  and  presentation  to  a  decision  operator. 

This  particular  problem  of  random  target  signals  will  not  be  pursued 
further. 

The  system  depicted  in  Fig.  4.10  is  completely  linear  and  produces 
the  best  peak  target  signal  power  to  root  mean  square  background  signal 
possible.  However,  a  complete  alarm  system  also  contains  a  threshold 
decision  filter  which  is  inherently  non-linear.  An  alarm  is  sounded 
if  a  signal  is  found  to  exceed  a  given  threshold  value  or  some  other 
condition  is  fulfilled.  There  is  usually  some  finite  probability  that 
the  background  alone  will  exceed  the  threshold,  and  it  is  just  the  peak 
background  signals  that  really  concern  us.  One  cannot,  therefore,  expect 
the  linear  system  which  maximizes  the  peak  target  signal  to  the  root -mean- 
square  background  signal  to  be  the  best  for  discrimination  purposes. 
Considerations  of  other  decision  criteria  and  discrimination  schemes 
lead  to  similar  conclusions.  An  example  of  another  criterion  is  developed 
in  the  following  paragraph. 

The  total  system  transfer  function  of  a  matched  linear  filter 
is 

K(a^)  =  (u>)  K2  (4.120) 

where  we  have  seen  that  Kj^(cO)  is  related  to  the  reciprocal  of  the 
input  background  spectral  density  and  K2(t*J)  is  related  to  the  input 
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signal  Fourier  transform.  Suppose  then,  that  a  detection  system  having 
such  a  transfer  function  is  required  to  distinguish  missile  targets  from 
severe  cloud  backgrounds  using  a  signal  gradient  criterion.  The  target 
is  associated  with  high  radiance  gradients  whereas  a  cloud  is  considered 
to  produce  much  smaller  gradients.  However,  the  transfer  function  of 
the  prewhitening  filter  raises  the  spectral  density  of  the  background  at 
high  frequencies,  making  discrimination  at  this  end  of  the  spectrum  more 
difficult.  The  advantage  of  the  point  source  resemblance  of  the  target 
is  at  least  partially  compromised. 

A  system  providing  the  maximum  ratio  of  target  signal  power 
to  root-mean-square  background  power  may  not  be  best  for  purposes  of 
discrimination  of  target  signal  presence  from  severe  backgrounds.  Indeed, 
we  have  considered  likely  cases  of  background  which  have  increased  proba¬ 
bility  density  above  the  value  6f  a  Gaussian  prediction  at  high  radiance 
levels.  Thus,  the  performance  of  the  optimum  linear  filter  could  be 
exceeded  by  a  system  specifically  designed  to  discriminate  on  the  basis 
of  some  characteristic  possessed  by  the  high  level  background  signals. 

Since  the  infrared  sensors  and  preamplifiers  are  inherent  parts 
of  infrared  detection  systems,  it  will  be  germane  to  consider  the  require¬ 
ments  placed  on  these  elements  by  the  nature  of  the  decision  operator. 

For  a  discrimination  scheme  which  requires  that  some  characteristic  of 
the  signal  be  preserved  the  minimum  mean-square  error  criterion  is  appli¬ 
cable.  For  systems  utilizing  correlation  as  discussed  in  the  previous 
section,  sensitivity  may  be  the  most  important  factor  which  suggests  the 


Actually,  the  requirement  that  the  probability  density  functions  of  target 
and  background  inputs  remain  essentially  unchanged  may  be  a  more  stringent 
requirement  than  high  sensitivity.  This  would  require  a  linear  system 
with  a  minimum  system  time  constant. 
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use  of  the  maximum  signal  to  noise  criterion.  Thus,  such  a  criterion  is 
by  no  means  inapplicable  to  our  problem  but  must  be  restricted  to  those 
subsystems  where  its  application  is  appropriate. 
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APPENDIX  4A 


WIENER  SPECTRA  ON  THE  SURFACE  OF  A  SPHERE 


This  appendix  is  concerned  with  the  description  of  the 
correlations  between  individual  random  variables  representing  the 
radiance  of  particular  resolvable  points.  The  usual  method  of 
description  by  autocorrelation  functions  or  Wiener  spectra  defined 
on  a  plane  is  very  arbitrary.  A  flat  plane  is  a  poor  representation 
of  the  earth  outside  a  very  limited  area  and  is  similarly  poor  for  an 
apparent  object  surface  of  a  particular  equipment.  The  arbitrary 
choice  of  a  finite  plane  approximation  allows  representation  of  a 
random  pattern  by  a  two-dimensional  Fourier  series.  This  repre¬ 
sentation  is  adequate  to  reconstruct  the  pattern  on  the  plane.  It 
does  not  allow  correct  matching  of  two  partially  non-overlapping 
patterns  and  the  plane  approximation  itself  can  only  be  used  over 
a  small  area  where  the  difference  between  spherical  and  plane 
trigonometry  is  comparable  to  a  desired  spatial  resolution. 

Consequently,  it  was  decided  to  investigate  the  correlation 
properties  of  radiance  patterns  on  a  sphere  representing  the  earth 
itself.  The  distortions  involved  in  viewing  these  patterns  would  then 
be  considered  to  be  part  of  the  information  processing  in  the  observing 
system. 

It  was  first  shown  that  the  radiance  patterns  observed  with 
any  arbitrary  resolution  could  be  represented  (in  the  sense  of  recon¬ 
struction)  by  a  series  of  Laplace  spherical  harmonics.  The  asymptotic 
properties  and  correlation-power  spectral  properties  of  these  series 
were  investigated.  It  was  found  that  the  Wiener-Khintchine  theorem  is 
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not  valid  in  spherical  geometry.  The  principal  virtue  of  Fourier 
representation  is  therefore  nullified. 

It  is  suggested  that  this  observation  has  serious  conse¬ 
quences  for  Wiener  spectral  analyses  of  background  in  general.  Two- 
dimensional  plane  spectra  are  often  used  with  the  explicit  or  implicit 
justification  that  a  plane  is  an  adequate  approximation  to  a  sphere 
over  a  small  region.  It  is  assumed  that  the  properties  of  the  plane 
Wi  '  .er  spectra  could  be  extended  to  spherical  spectra  at  the  expense 
only  of  labor  if  the  extra  accuracy  were  ever  required. 

Suppose  that  the  near  infrared  radiance  pattern  of  the  earth 
is  observed  with  a  finite  angular  resolution  OC  from  a  series  of 
satellites  at  altitude  h.  Then  the  linear  resolution  of  the  measure¬ 
ments  at  the  earth  surface  is  never  less  than  CCh.  The  radiance  is 
also  bounded  in  magnitude  to  a  value  at  least  less  than  a  specular 
reflection  of  the  sun.  Consequently,  the  observed  pattern  is  con¬ 
tinuous  and  absolutely  integrable.  It  may  be  represented  by  a  function 
g(8,0)  where  the  coordinates  9  and  0  are  spherical  angular  coordinates. 
Absolute  integrability  implies  that  the  function  g  will  have  a  Fourier 
transform.  The  assumption  of  finite  resolution  implies  that  the 
transform  of  the  observable  process  will  be  band  limited. 

It  would  also  be  possible  to  have  a  sufficient  number  of 
satellites  spaced  around  the  earth  so  that  the  entire  surface  would 
be  observed  simultaneously.  In  this  case,  the  complete  function 
g(9>0)  must  obey  the  periodicity  constraints 

g(e.0)  =  8(6,  0  +  2  >r) 

8(6,0)  =  g(6  +  TT,  0) 

Because  of  these  constraints  and  band  limited  Fourier  trans¬ 
form  ability,  the  function  g(0,  0)  can  be  represented  in  all  the  detail 
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that  it  can  be  observed  by  a  Fourier  series  with  a  finite  number 
of  terms.  The  number  of  terms  will  be  of  the  order  of 
2  (  477'R^/(a:h)2  }  .  A  particular  form  of  Fourier  series  which 
can  be  used  to  represent  a  function  in  spherical  coordinates  is  the 
series  of  Laplace  spherical  harmonics  which  consist  of  products  between 
a  series  of  complex  exponentials  in  h0  and  associated  Legendre  poly¬ 
nomials  in  the  n,  h  and  0.  In  this  series  the  harmonic  of  order  h,n 
has  h  periods  around  the  0  =  'n'/2  (equator)  circumference  and  n 
periods  around  a  0  =  0  (meridian)  circumference.  The  series  expansion 
is  set  out  explicitly  in  the  analysis  section  below. 

Conceptually  a  class  of  background  radiance  processes  is 
selected  as  generated  by  a  weather  condition,  a  point  of  observation  and 
a  sun  location.  A  number  of  sample  functions  of  this  process  are 
observed  on  different  occasions  from  an  array  of  satellites.  Each 
sample  function  is  broken  down  into  its  representation  as  a  series 
of  Laplace  spherical  harmonics.  Thfc  coefficients  c^^^  of  this  series 
together  contain  enough  information  to  reconstruct  each  sample  function. 
Application  of  this  technique  to  a  number  of  sample  functions  will 
yield  a  set  of  values  for  each  The  distribution  of  this  set  of 

values  is  an  estimate  of  the  probability  distribution  of  the  coefficient. 

4A. 2  SPHERICAL  HARMONICS 

Any  particular  radiance  pattern  on  the  sphere  can  be  repre¬ 
sented  by  the  series 


g(e.0) 


■I 

n^ 


no 


P  (cos0) 
n 


(Snh  cos  h  0  +  b^^  sin  h  0)  P^  (cos  0) 


nh 
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where  the  a,  and  b  ,  are  a  series  of  Fourier  coefficients,  the 
nh  nh 

h 

are  associated  Legendre  polynomials  of  the  first  kind.  The 

series  can  be  put  in  a  more  compact  form  as 

OO  n 


g  (9,0)  =  Y.  Z 


nK)  h-n 


because 


a  cos  h  0  H-  b  sin  h  0  =  ^  a  ,  (e^^^  +  e 

nh  nh  /  nh 

1^0  2 


,  nh  .  ih0  -ih0v 

+  IT  ("  -  "  ) 


11 

=  1: 


C  U  ® 

nh 


where  c 


s  u  “  u 
nh  nh 


for  h  >  0 


c  ,  =  a  ^  for  h  =  0 

nh  nO 


^  u  +  u 

nh  nh  ^  ^  ^ 

c  ,  - r -  for  h  <.  0 

nh  2 
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* 

^nh 


'n-li 


The  ct)e  ri'ic  i  ents  are  found  f  roin  the  iuLej^^ra  I  s 

.'tr 

( 

2n  +  I 

-yr  0 


no  “■  47r 


r  r 


sin  0  d0  g  (0,  0)  (cos  0) 


/V  -TT 

"nh  °  "4^  ^  0  d0  g  (0,0)  P_  '’  (cos  0)  cos  h0 


/4r 


2n+l  (n-li)  ! 


nil  4/r  (n+li).' 


2n+l  (n-h) ! 
^nh  ^  ^nh  (n+h)  ! 


;  /  d0  |sln0ci0  g(0,0)  I’  (cos  0)  sin  h0 


d0  /  sin  0d0  g  (9,0)  P  ^(cos  0)  e 

*  n 

-TT  0 


4A.3  CHANtilNC  COORDINATKS 

In  orde  r  to  i  lives  t  ij^a  te  Llie  description  of  correlation  by 
these  series  it  will  be  necessary  to  transf(3rm  coordinates  0,0  to 
others  0^,  0^  through  wliat  ainouiils  to  a  rotation  of  the  sphere.  The 
angle  of  roiat ion  is  jL'iven  by 

cos  cx  =  cos  0  cos  0^  +  sin  0  sin  0^  cos  (0-0 S 

l>ecause  the  scalar  product  oi  the  vectors  (  1,0,0)  and  (  1,0^, 0^ 

A  1  1  1 

A. A  =  XX  +  yy  +  xz 
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=  cos  0  cos  0^  sin  0  sin  0^  +  sin  0  sin  0^  sin  0  sin  0^  +  cos  0  cos  0^ 
=  cos  0  cos  0^  +  sin  0  sin  0^  (cos  0  cos  0^  +  sin  0  sin  0^) 


cos  0  cos  0^  +  sin  0  sin  0^  cos  (0-0^) 


But  A. A  =  cos  oc  ^ 

It  can  then  be  shown  that  P  ^  (cos  0^)  e^^^  can  be  expressed  in  terms 

n 


of  0  and  0  and  the  rotation  required  to  carry  0,0  into  0^,0^  by 

n 

p;  (cos  .s  .  z  c  p;  (CO.  « 


r=-n 


where 

ohr 


2n  (n+r)  .* 


•i(h+r)/>  ^-i(r-h)<r  (cosr)^'"’'  (sin 


,n+r 


/  2^.  n-h  .  2  -vn+h 

2  r-^> 


d (cos  Y' ) 

and  the  and  'ip  are  defined  in  terms  of  the  direction  cosines 
of  the  axis  of  rotation  from  0,0  to  0^  0^  and  the  magnitude  of  the 
rotation  as  follows 


Direction  cosines  of  axis  of  rotation  are 


w. r .  t. 


-1 


X  ^  sin^  sinT  ^  =  tan 


"  y  ^  cos«r-'sin  T  /^  =  tan'^  (^3/^4) 


qs/ti)  sin^  cos  sin  ^ 
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Magnitude  of  rotation  OC  is  related  to  q. 


cos  ^  **  *14  “  cos  p  cos 


1/2^2 

cos  ^3  ^4 

2 
V 


and 


.  ^ 


sin 


=  v  >  ^  =Vl^  ■*■ 


2  ^  2 

+  q3  ;  V 


2  .  rr". — T'.  r 


■  ’1^  \f 


2  '^‘13  \ 


specifically  if  v  =  1  ,  U)  =  sin  CCl2, 

The  direction  cosines  of  the  axis  of  rotation  may  be  determined  as 
follows : 

Let  an  initial  vector  A  be  (1,0,0) 
a  lagged  vector  B  be  (1,0^,  0^) 

A  =  cos  0  sin  0,  sin  0  sin  0,  cos  0 


B  =  cos  0  sin  0  ,  sin  0  sin  0  ,  cos  0 


AxB=(AB  -  AB),(AB  -  AB))(AB  -  AB) 
yz  zy'^zx  xz  xy  yx 


The  axis  of  rotation  unit  vector  is  (A  x  B).  (i,  j,  k)/sino^ 


(qi,q2.q3)/sin^  =  ^ 


(sin  0  sin  0  cos  0^  -  cos  0  sin  0^  sin  0^)  , 


(cos  0  cos  0^  sin  0^  -  cos  0  sin  0  cos  0^), 

(cos  0  sin  0  sin  0^  sin  0^  -  cos  0^  sin  0^  sin  0  sin  0)^/sin^ 
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Solving  O^,  f  ,X' 


tanO^«  - 


_  sin  0  sin  9  cos  -  cos  Q  sin  0^  sin  9^ 
^2  cos  0  cos  0^  sin  0^  ~  cos  0  sin  0  cos  0^ 


\P  =  —  =  (cos  0  sin  0  sin  0^  sin  0^  -  cos  0^  sin  0^  sin  0  sin  0) 

I 


.1  ... 


sinO^/2 _ 

sinOC  cos  ^/2 

2  1  112 
sin  »  (sin  0  sin  0  cos  0  -  cos  0  sin  0  sin  0  )  + 

11  12 
(cos  0  cos  0  sin  0  -  cos  0  sin  0  cos  0  ) 

cos^  It  =  (cos  0  sin  0  sin  0^  sin  0^  -  cos  0^  sin  0^  sin  0  sin  0)^  +  cos^^ 

=  (cos  0  sin  0  sin  0^  sin  0^  -  cos  0^  sin  0^  sin  0  sin  0)^ 

11  11  12 
(cos  0  sin  0  cos  0  sin  0  +  sin  0  sin  0  sin  0  sin  0  +  cos  0cos0  )  -1 

- 2-  - - - - 

22221  1  11 
sin  T'  =  sin  0  sin  0  cos  0  -  2s  in  0s  in  0  cos0  cos0  sin0  sin  0  + 

2  21  21  2  2  21 

cos  0  sin  0  sin  0  +cos  0  sin  0  cos  0 

1  1  12  2121 
2  cos  0  sin  0  sin  0  cos  0  cos  0  cos  0  -fcos  0  cos  0  sin  0 


=  sin  0  cos  0 

-  2  sin0  cos0  (sin0  cos0^  sin0^  sin0^  +  cos0  sin0^  cos0^  cos0^) 

2  2  1 
+  cos  0  sin  0 
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221  22l  .11 

sin  0  cos  0  +  cos  0  sin  0  -  2  sin0  cos0  sin0  cos0 

(sin0  sin0^  +  cos0  cos0^) 

sin^0  cos^0^  H-  cos^0  sin^0^  -  2  sin0  cos0  sin0^  cos0^  cos  (0  -  0^) 


Note  if  0  =  0 


I  1  12 

sin  ^  =  (sin0  cosG  -  cos0  sin0  ) 


=  sin^  (0  -  0^)  . 

r  =  +  (e  -  +  2  ^ 


tancr=-  =  tan  0  O' -  0  +  2 

COS0  (cosG  sin0  -  sinG  cosG  ) 

tan^  sin0  (sin0  sinG^  -  sin0  sin0^)  _  ^  0  +  nTT 

2  cos^  CC 


Also  note  that  if  0^  =  0 
,1 


0  =  0  +7| 

2  2^  2  2  2  2  ^  1 
in  Jf'  =  sin  0  /cos  0  cos  ^  +  sin  0  sin  ^  -  2  cosG  sinG  sin^J 

2^22A  22  ^  ^ 

+  cos  0  ^  sin  0  cos  £  +  cos  0  sin  ^  +  2  cos0  sin0  cos^ 

-  2  sinG  COS0  cos0  sin0  cos^^  -  sin^G  ^  sin^H- 
0  cos  ^  “■  sinG  sin^^J  cosj^ 


cos 


2^  1  .  2  y, 

cos  ^  *  1  -  sin  Z 
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4A,4  AUTOCORRELATION  AND  POWER  SPECTRA 

For  a  one-dimensional  Fourier  series  representation  on  an 
interval  (-T,T) 


.00 


i  2  ?7nx 

Gn  e 


r 

3  * 

n 

Jn 


n  =  -00 

T  -I'TT'nx 


G  *  I  g(x)  e 
‘'-T 


dx 


the  two  quantities 


R(^  )  -  2t 

T-^00 


Lim 


I 


g(x)  g(x  +  dx 


and 


Lim 


G  G 
n  n 

2T 


may  be  defined.  They  are  called  the  autocorrelation  function  and  power 
spectral  density  of  the  random  process  g(x).  They  have  far  reaching 
statistical  implication  if  the  process  is  stationary.  Substituting  the 
series  for  g(x)  into  the  integral  for  R(^ ) 


Lim 


i2^  nx 

e  2T  dx 


Interchanging  the  orders  of  limiting  processes  (which  is  permissible  if  g(x)  is 
absolutely  integrable)and  changing  variables  to  y  =  x  +  ^ 
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“(f)  ’X. 


ilTt 

2T 


■H 


Lim  -2^  g(y)  e 

T400  /  . 

-T  +5 


dy 


since  the  process  g(x)  is  stationary  this  is 


+12^  ni 


which  also  iraplies 


W  =  Lim 

n  T->«» 


>(f) 


-  i2  TT 

- 2T - 


-T 


This  Fourier  transform  relationship  between  R(^  )  and  is  called  the 
Wiener-Khintchine  theorem.  It  is  the  purpose  of  this  analysis  to  see 
whether  a  similar  relationship  holds  for  the  Laplace  spherical  harmonics. 

4A.5  CORRELATION  ON  A  GREAT  CIRCLE 

Consider  first  two  special  cases 

(a)  g(e,  0)  *  g(0) 

and 


ih0 


(b)  g(e,  0)  =  g(0) 

oo 

In  case  (a)  we  have  g(0j  0)  =  ^  C^ 

h= 

and  clearly  x  has  been  replaced  by  0  with  no  other  change  but  that  the 
range  of  0  is  restricted  to  lie  in  (-  it, IT).  The  limiting  process  specified 
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above  is  meaningless  but  the  other  identities  carry  through  formally. 
Case  (b)  may  be  written  as 


g(6>0)  =  g(cos  6)  =  ^ 

n*0 

2n+l 


(cos  0) 


^n=  2 

-tT 


1 


g(cos  0)  (cos  0)  sin  0  d0 


Now  define  R( 


»  I  sin  0  d0  g(cos  0)  g  ^ cos  (e 


and  substitute  for  g/cos(0  “f  )( 

rlT  ^  ^  ^ 

R(^)  = 


sin  6  d0 


oo  ^ 

g(cos  0)  ^  ^ cos  (9-0  I 


Pn  003(9-^  )]  =  (SS)’'  C  ®>  C 

cos  m  (0-01) 

taking  0=0^  and  m  =  0  which  amounts  to  assuming  that  g(0,0)  is  constant 
around  lines  of  latitude 


Pn  •  0]-^.  (cos  0)  P^  (cos^) 


R( 


r 

0  ~  ^n  I  0  d0  g  (cos  0)  P^  (cos  0) 


=  ^  P^  (cos^) 


which  is  a  form  of  the  required  theorem. 


(cos 
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4A.6  CORRELATION  IN  SPHERICAL  HARMONICS 


In  the  general  case  define 

V  .If 

R(£,7[)  =  j  d0  J  Sin0  d0  g(e,0)  g(e\0S 

-'If  '^0 


with  the  intention  of  showing  that  this  is  equivalent  to 

n 


CO 

)  =  constant  ^ 


c  ^  gih(0-0  )  p  h  (0-0^) 

nh  n  ^  ^ 


n=0  h=-n 


=  constant 


E  E  f  > 


n=0  h-n 


R  (^^)  =  I  d0  j  sin0  d0  g(0,0)  ^  ^  Cnhe^^‘**  P^^(cos0^) 

I  n=0  h=-n 

fd0  sine  do  g(0,0)][]  Ec„i^  Z  ^T^i'  ®2n 


n  h  r=-n 


-E  E 


n=0  h=“n  A 


n 

V*  (n-  r  ) !  r  hr  f  , 

h  J  j  “ 


in0d0g(0,0)  P  (co80) 
n 


The  algebra  required  to  perform  this  integration  and  compare  with  the 
desired  result  is  very  tedious  in  general.  It  is  therefore  worthwhile 
to  test  a  special  case  in  order  that  the  labor  may  not  be  wasted. 
Therefore,  try  a  simple  case  where  n  =  2  and  0=0^,  (Note  that  0=0^ 
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does  not  imply  a  great  circle  correlation  necessarily). 
2 

g  (6.0)  =  Yj  P2^  (cos  0) 

h=-2 

all  other  coefficients  are  zero. 

Also  sin^^  »  sin^(0-9^)  =  sin  ^  ^ 
tan  CJ^  w  tan  0 
tan  =0 

Taking  the  principal  solutions 


hr 


(2+r  )  ; 


•  /  uNrf  h+  r  r -h  ,2+  r 

e  ^  ^  (cosT)  (sin  7"  )  - 


d(cos  7-) 


,  2  .,2-h 

r„2+r 


.  .  2^.2+h 
(sin  X) 


SO  that 

[-‘■•■rj 

^  r 

sinG  d0  g(9,0) 

>'  hr  ^ 

^  2n 

2lf 

(n+h/: 

“  ^-i(r-h)0 

2n+l 

(n-h)  .' 

and 

2 

R(£>7))  =  4-rf  ^ 

*  n=2  h= 

'-2  , 

in  this  case 


s  hr 

(2+h)!(2-r  )■'  _4 


(2-h) 


.2 


-i(r-h)0 


C 

2h 


and  it  is  to  be  seen  whether  this  is  equal  to 


We  shall  proceed  to  tabulate  and  compute  these  summations  taking  the 
simpler  second  one  first. 
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TABLE  1 


TABULATION  OF 


h~2 


^  X3fa 

h^2  2^2 


h 

2 


<'2  2 
Tl+iT  <y 

dy 


Term  Coefficient  (l-y  )  dv 


h  2+h 

2  — fevTCy  -1)  Value  of 

Derivative 


Total 

Coefficient 


C22  1/8 


(i-y^) 


d  /  2  ,.2 

—4 (y  -1) 

dy 


24 


3(l-y^) 


■=21 


/I  2A  d^  ,  2  ^.2 
(l-y  y  — 5  (y  -1) 

dy 


24y 


3y  \^-y^ 


Cjo  1/8 


„  2.0  d‘  .  1  -.2 

(l-y  )  — j  (y  -1) 

dy 


4(3y^-l)  j(3y^-l) 


S-1 


„  2.-^  d  ,  2  ,.2 

(l-y  )  ^  (y  -1) 


4y(y  -1)  -%yYl-y 


^2-2 


2.-1  .  2  -.2 

(l-y  )  (y  -1) 


/I  2.2  .  1  ,,  2. 

(l-y  )  +3  (l-y  ) 


Consequently  the  entire  series  should  be 

3(i-yW^2  +  3y  '^^^21  i<3y^-i)  C20  -  fy  +  |(i-y^)  cl_^ 

-  3  sln^£C22  +  3cos£sln^C2i  +  |(3co8^£-  l)C2o=-|cos£ sln^C^,^  +|sln^£c 
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TABLE  2 


f  27r 

r?.2  ^ 

(2+h).'«-r);  ' 

(2-h) 

Li!  (co»£) 

^-i(r-h)^ 

-  l/» 

S2:^L 

<2-l>).'^ 

.r-h 

2  .2+r 

Putting  a 

■  co«*^ 

h  t 

h  coaff 

r  coaff 

h»r  r-h 

*  ^  a-i)  ^ 

Darivativa 

Valua  of 
Darivativa 

Product 

-2 

2 '10^15 

4.' 

<1-)° 

da® 

.*  <.-«* 

2‘lf'li^ 

-2  -1 

3.' 

da 

4.’ 

-2  0 

2'70u:5 

21/21 

«■*  (i-.)‘ 

ii 

da^ 

12a^ 

-2  1 

27?74.‘5 

1/3.' 

da^ 

24c 

-2  2 

2V'/a:5 

1/4.' 

d^ 

24 

Total  tam  in  ^oa^  ^  ^  ^ 

+  ?oo.£.lii’^ 

:  *54 •*““£) 

-1  -2 

2'V'  1 

T  31 

4.' 

(I-.)-* 

i_ 

da® 

a^(a-l) 

.  ^<4 

-1  -1 

Ilf  1 

T  IT 

3.’ 

(1-.)* 

da 

a^(4aO) 

V?^.(4.-» 

•1  0 

1 

.■*  (1-.)' 

i_ 

da^ 

6a(28.1> 

-1  1 

2'Tf  1 

T  IT 

1/3; 

••  (1-)' 

d.’ 

6(4a.l) 

0-.)(4.-l) 

-1  2 

211'  1 

5  IT 

1/4: 

dm* 

24 

ao  that  tha  total  tan  in 

ia 

^  2- 

«  4  coa' 

•iL 

coa^  ^ (4  coa^^ 

^  -  3)  +  coa£aln^(2 

-1) 

+  i  ain 

£  (4  coa^, 

£-1)  +?«•£ 

•‘“’£j 

-ULc* 

5  ''2-: 

1  ♦ 

-  2  coa^^  ( 

il«£+  coa^^(^ 

ain  £-3)  +  coa 

£ita^<i  .te* 

£-‘)  -5“*£] 
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It  is  now  clear  that  the  terms  in  C2.2  ^2-1 

on  Table  2  have  different  values  than  those  obtained  on  Table  1.  For 
example,  when  ^  — >  0  the  term  in  ^2^-^  goes  to  zero  on  Table  1  but 
to  277^/5  on  Table  2. 

The  values  of  the  describe  the  radiance  pattern.  They 
can  be  chosen  arbitrarily  from  any  set  which  keeps  the  pattern  positive, 
bounded  and  continuous.  Consequently,  the  discrepancies  in  the  function 
of  cos  ^  are  not  to  be  corrected  by  relations  between  the 

Thus  the  Wiener-Khintchine  theorem  does  not  hold  for  the 

particular  spherical  harmonic  series 
2 

g(e,0)  =  (cos  0) 

h=-2 

when  the  lag  0-0^  is  arbitrary  but  0-0^  is  zero.  Consequently,  the 
theorem  does  not  hold  in  general  for  spherical  harmonics. 
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Attn:  Dr.  C.  Cummlng 

1 

Aerospace  Corporation 

2400  El  Segundo  Boulevard 
£1  Segundo,  California 

Attn:  Dr.  C.  Shervln 

H.  Wesseley 

I.  Spiro 

3 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Attn :  Dr .  S .  Passman 

Dr.  D.  Delrmendjlan 

2 

Denver  Research  Institute 

University  of  Denver 

Denver,  Colorado 

Attn :  D .  Murcray 

D-3 

1 

Commander 

Air  Force  Cambridge  Research  Center 
Geophysics  Research  Directorate 
L.  G.  Hanscom  Field 
Bedford y  Massachusetts 
Attn:  CRZC 

Commander 

Aeronautical  Systems  Division 
Air  Force  Systems  Command 
Wright-Patterson  AFB,  Ohio 
Attn:  ASRNS 

Commander 

Air  Force  Space  Systems  Division 
Air  Force  Unit  Post  Office 
Los  Angeles  45,  California 
Attn:  Major  J.  Wallace 

Colonel  Norman  (461  Office) 


